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ABSTRACT 



This thesis examines spectrometric oil analysis data 
from two sources in an attempt to formulate a statistical 
model which will be useful in monitoring aircraft engines 
in the Naval Oil Analysis Program. Initially , experimental 
data, gathered for an Air Force study, is used to determine 
if the measurement error inherent in the monitoring 
procedure is normally distributed and if correlations exist 
between measurements for different wear metals. Based on 
the results of this investigation, a study is made of 
operational data from Wright reciprocating engines of the 
R1820-82 model type. This investigation leads to the 
conclusion that a multivariate regression model is useful 
in estimating the parameters of the distribution of analyses 
from properly operating engines of this type. A procedure 
is then suggested which would employ the readings from past 
oil analyses from a particular engine to determine its 
present condition. 
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I. 



INTRODUCTION 



For several years the technique of using the spectro- 
metric analysis of oil samples as an aid in determining the 
condition of diesel engines has been employed successfully 
by major railroads and various other users of large diesel 
equipment. In 1956, a trial program was begun at the Naval 
Air Rework Facility in Pensacola to determine if this 
method could also be used to monitor aircraft engines. 

Since that time the program has proved successful and has 
evolved into the Naval Oil Analysis Program (NOAP) . It is 
planned that this program will eventually include all Navy 
fluid lubricated mechanical systems. A more detailed 
history of the program is contained in Refs. 1 and 2. 

Since this thesis is concerned with an investigation of 
data collected at the Pensacola laboratory, the following 
descriptions of the operation will be limited to the pro- 
cedures used there. Reciprocating aircraft engines are 
sampled approximately every 30 hours. The sample is taken 
after the aircraft has returned from a flight and before 
the oil has become cold. It is immediately sent to the 
laboratory by air mail and is analyzed on the day received. 
The analysis is accomplished by a spectrometer using the 
rotating graphite electrode technique. Measurements of the 
parts per million (ppm) content of ten metallic elements, 
which might be indicative of engine wear, are made 
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simultaneously. Of these ten, aluminum, chromium, iron, 
silver, copper, magnesium and nickel are those which are 
relevent to engines of the model considered in this report. 
The ppm readings are automatically recorded on a punched 
card which also contains various other hand-entered data 
identifying the sample. 

Once the data has been recorded, it is used to aid in 
determining what the operating condition of the engine 
might be. Presumably, if the engine is operating properly, 
the amount of metallic contamination in the circulating oil 
should be within certain normal limits. In addition, it is 
felt that the amount of contamination added to the oil 
since it was last sampled should be within specific limits 
if the engine is in good working condition. If, however, 
the engine is discrepant and excessive wear is present, 
this will presumably cause an abnormal addition of metallic 
contaminates to the circulating oil. 

Thus, when a sample has been analyzed and the results 
recorded, both the magnitudes of the present readings and 
the changes in the readings since the last sampling are 
compared with threshold limits which have been developed 
for each engine type and each metallic element relevent to 
that engine. If the results fall outside the prescribed 
limits, some action is generally taken by the laboratory. 
Usually another sample is requested and the previous 
results are verified. If the abnormality persists, either 
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the aircraft is grounded for maintenance or future samples 
are taken more frequently than the usual 30 hour interval. 

At present, these threshold limits are subjectively set 
and vary only from element to element and among engine 
model types. They are based on the past history of the 
aircraft model which includes the data supplied by the 
engine manufacturer before the model is placed in service 
and experience accumulated once the model is in use. The 
limits are not used as sharp boundaries for classifying 
engines as normal or discrepant but merely as indicators 
upon which a subjective decision as to the action to be 
taken can be based. 

This report examines two sets of data from the Pensa- 
cola laboratory with the intention of determining the 
propriety of three assumptions implicit in this classifi- 
cation procedure. Since the same threshold limits are 
used for all engines of a particular model, it is assumed 
that all normally operating engines of the same type can be 
expected to have the same amounts of metallic contamination 
in their oil systems. In addition, since threshold limits 
are constant for a given element and model type, variations 
in other factors, such as the operating hours since the 
last oil change, must be ignored or subjectively introduced 
into the classification procedure. Finally, since thresh- 
old limits are set for each element independent of the 
limits for other elements, readings for different elements 
are assumed to be uncorrelated. 
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Once these assumptions are verified or rejected, a 
statistical model is formulated to aid in establishing a 
more objective classification criterion. 
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II. ERRORS INHERENT IN THE MONITORING PROCEDURE 



Since the intention of NOAP is to make inferences 
about the condition of aircraft engines , based on the 
amount of wear metal contamination in the engine's oil 
system, it is extremely important that the amount of con- 
tamination recorded at the laboratory accurately reflect 
the actual amount present in the engine. For the purposes 
of this report, measurement error will be defined as the 
difference between the parts per million content of an 
element recorded as present in a particular engine at a 
point in time and the actual content at that time. In NOAP 
there are a variety of potential sources of error, all of 
which can contribute to the net measurement error defined 
above . 

A. ERRORS IN SAMPLING 

As was mentioned earlier, oil samples are taken from 
reciprocating engines normally every 30 flying hours and 
while the oil is still hot. This sampling is accomplished 
with a special sampling kit consisting of a sampling tube 
and a sampling bottle. The tube is inserted into the oil 
reservoir, and when it has filled the top end is stopped 
with the operator's finger. The contents are then trans- 
ferred to the bottle, which is immediately forwarded to 
the laboratory for analysis. When the sample is analyzed, 
a small portion of the oil in the bottle is used in the 
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analysis [Refs. 1 and 2], Thus, an extremely small amount 
of oil is used to determine the extent of contamination in 
the engine's entire oil system. Any lack of homogeneity 
in the engine's oil reservoir will result in a non- 
representative sample. Further, any contamination added 
to the sample through a lack of cleanliness of the sampling 
tube and bottle or through handling at the laboratory will 
contribute to the measurement error. 

B. ERRORS IN RECORDING 

At the time the sample is taken, certain data including 
the date, the operating hours since the last oil change, 
the hours since the last overhaul of the engine, the engine 
serial number and the model number are recorded and mailed 
to the laboratory with the sample. Various portions of 
this data are transferred from other records. At the 
laboratory the data are entered by hand on the permanent 
record cards maintained there [Ref. 1] . This entire 
sequence of recording and transferring data from one record 
to another can result in mistakes. 

Unfortunately, as with the errors in sampling, there is 
no data available at the present time that can be used to 
measure this error. 

C. ERRORS IN ANALYSIS 

When an oil sample is received at the Pensacola labora- 
tory, it is analyzed using a direct reading spectrometer 
with spark excitation, stationary and rotating disc 
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electrodes. The sample bottle cap is filled with oil and 
placed in the spark stand. The gap between the two 
electrodes is set and the disc electrode begins to rotate 
at 30 rpm. As the electrode rotates, a thin film of oil 
is forced to the area under the fixed electrode. A high 
energy spark is then fired across the gap and the film of 
oil is burned for 25 seconds. The light from the burning 
oil is separated so that its intensity at the wave lengths, 
produced by the elements to be analyzed, can be compared 
with built-in standards. The average intensity over the 
burning period is then measured for each element simul- 
taneously and converted into parts per million. These 
readings are automatically recorded on the engine history 
card [Ref. 2] . 

If it is assumed that there were no errors in the 
sampling or recording and thus, that the oil used in the 
analysis is representative of the oil in the engine, any 
difference between the true content of contamination in 
the engine and that recorded after the analysis can be 
attributed to an analysis error. An experiment designed 
to measure the effects of this type of error has been 
conducted and the results of an analysis of this data are 
presented in the next section. 
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III. DISCUSSION OF THE ANALYSIS ERROR 



Although the error due to the spectrometric analysis 
of the oil is not the only possible source of error, it 
certainly is a major contributor to the over-all measure- 
ment error defined earlier. For this reason, an examination 
of data accumulated for a study conducted by the Air Force 
[Ref. 3] was performed and the results are discussed in 
this section. 

A . DATA 

In 1967, the Pensacola laboratory participated in an 
experiment conducted by the Air Force. During a 30 day 
period the laboratory received 100 oil samples. These 
were to be analyzed in the normal manner and the results 
reported. Although the laboratory was not aware of it, 
these 100 samples consisted of ten samples each repeated 
ten times. Thus, the laboratory actually repeated the 
analysis of ten different samples ten times. The results 
of the analyses, for the seven elements of interest in 
this report, are included in the Computer Output section, 
where the readings on like samples are in groups numbered 
from one to ten. Missing data accounts for some groups 
having less than ten readings. 

For each of the ten groups of repetitious analyses, the 

sample mean and standard deviation were calculated for each 

th 

element. If, for example, is the i — reading for 
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aluminum in a group of size n, then the sample mean, X, 
and standard deviation, S, for aluminum are 



n 

X = ( Z X. )/n 
i=l 1 



and 



n _ 2 

S = [ Z (X. - X) / (n-1) ] 
i=l 1 



( 1 ) 



respectively. The results of these calculations are also 
presented in the Computer Output section for each element 
and each group of samples. 

In addition, for each sample group, an estimated corre- 
lation matrix, R, was calculated, where if r. . is the 

— 

element in the i— row and j— column of R, then 



n 

Z 

k=l 



Z (X. , - X. ) (X . , - X . ) 

i,k i : ,k y 



ID 



n 

[ 2 (X 
k=l 



- 2 n - 2 Z 
. ,-X.)^ Z (X., - X . ) ^] 

i v-i D , k D 



( 2 ) 



where X. , is the k — reading for the i — element, and X. 

1 f K 1 

is the element's sample mean. These correlation matrices 
are included in the Computer Output section. 

These preliminary computations provided statistics 
which were used to test certain hypotheses concerning the 
probability distribution of the analysis error. 
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B. TEST FOR NORMALITY 

Since the overall measurement error is the net effect 
of errors arising from a variety of sources , the Central 
Limit Theorem of Probability Theory [Ref. 4] provides good 
justification for making an assumption of normality in the 
distribution of this error. Thus, any additional evidence, 
which tends to point to the normality of one of the con- 
tributing sources of error, will serve to strengthen the 
overall assumption. The experiment conducted by the Air 
Force provided data which was used to test for normality 
in the distribution of the analysis error. 

If is defined as a seven-component vector of sample 
readings arising from sample group k = 1,2,..., 10, then 



4 = Hk + 



where is a seven-component vector of the true metallic 

content of seven elements in the sample associated with the 
th 

k — group of readings, and e^ is the seven- component random 
analysis error vector. Thus, if it is assumed that e^ is 
a multivariate normal random variable with zero mean vector 
and unknown covariance matrix , then is a multivariate 
normal random variable with mean vector and covariance 
matrix , denoted N (jj^. , Z^) . 

If this assumption is correct, it is possible to perform 
a transformation of the form, 



Z k = P k (X k - V' 
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which will produce a multivariate normal random variable Z 



which has mean vector zero and covariance matrix I, the 
identity matrix. For details of this transformation see 
Appendix A. 

For each of the ten groups of sample readings the mean 
vector was estimated using 

n k 

= . E *kj /n k (3) 

:=i J 



where X^j is the j — vector of sample readings from the k — 
sample group, and n^. is the number of readings in that 
group. In addition, the covariance matrices, , were 
estimated using the unbiased estimator 



E. = (S . S . r . . ). 
— k 1 3 13 k 



(4) 



th 

where S. is the estimated standard deviation for the i— 

1 

element, computed as in equation ( 1 ) , and r^ is as defined 
by equation (2). For each of the ten groups the non- 
singular matrix, P^, was found and the transformation. 



Z. . = P. (X, . - X. ) , 

— k,i — k — k,i — k' ' 



performed on each vector of readings, X, ., in the k — 

k , 1 

group, k=l,2,...,10. In this way vectors Z, . were 

K f 1 

produced, the components of which are stochastically inde- 
pendent and are distributed according to N(0,1) if the 
hypothesis is true. 
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All readings from the ten groups were then pooled to 
produce a sample of 651 deviates , assumed to be univariate 
normal. The Kolmogorov- Smirnov goodness of fit test was 
applied to this sample and the resulting test statistic of 
.0215 was not significant even at the .20-level. Thus, the 
hypothesis of normality in the distribution of the analysis 
error was accepted. For details of the test see Appendix A. 

C. TEST FOR EQUALITY OF COVARIANCE MATRICES 
Let again be defined as 



as in the previous section, where now e^ is assumed to be 

a multivariate normal random vector. In addition, let the 

/\ 

estimators and Z_^ be defined by equations (3) and (4) 

respectively. Then, the Air Force data can be used to 

/\ 

produce ten estimated covariance matrices , each asso- 
ciated with a different sample group and thus, a different 
true content vector y^. If the true covariance matrices, 
are independent of the vector , and thus constant 
for all k = 1,2,..., 10, the ten estimated matrices could 
be pooled to obtain an over-all estimate of Z_. This 
hypothesis of the equality of the ten covariance matrices 
was tested, and the results led to the rejection of the 
hypothesis at the .10 level of significance. The details 
of the test used and the results obtained are included in 
Appendix A. 
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In an attempt to account for the apparent variability 
among covariance matrices from different sample groups, a 
regression model was formulated. It had been suggested by 
Baird-Atomic Inc. , the manufacturer of the spectrometer 
used at the Pensacola laboratory, that the variability in 
readings for a given element is dependent on the true 
content of the element. Specifically, the relationship 
is assumed to be 

2 . 2 
a = a + by 

2 

where a is the variance of repeated analyses of the same 
sample for a given element, y is the true parts per million 
content of the element, and a and b are constants [Ref. 5]. 
The model, 

S 2 = a + bX 2 + e, 

was used to examine the propriety of this relationship for 

2 

each of the seven elements under consideration. Here S 

is the square of the standard deviation estimate defined by 

_2 

equation (1), X is the square of the sample mean, and e is 
a random variable. For each element, a and b were estimated 
using least-squares techniques and, under the assumption 
that variations about the regression line are normally 
distributed, the hypothesis, b = 0, was tested [Refs. 6 
and 7] . A t-test was used and the significance level set 
at .10. Of the seven slopes tested, those associated with 
aluminum, iron, copper and magnesium were significantly 
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non- zero. It should be mentioned that the true content of 



the other elements did not vary much among the ten sample 
groups. The results of the least-squares estimation, 
together with the numerical results of the t- tests, are 
included in the Computer Output section. 

D. TEST FOR INDEPENDENCE AMONG ELEMENTS 

Because of the apparent dependence of the variance of 
repeated readings upon the true content of an element, it 
was felt that covariances between elements might depend 
upon the content of the elements concerned. If this were 
the case, then, for example, a particularly high reading 
of one element might be "explained" by a corresponding 
low reading of another element. Under the assumption of 
normality, the hypothesis of independence is equivalent to 
the hypothesis of zero correlation. This hypothesis was 
tested, using each of the ten correlation matrices R, 
defined by equation (2) . Of the ten tests conducted, only 
three of them were not significant at the .10 level. The 
numerical results and the test used are given in Appendix A. 
It can be noted from the correlation matrices in the 
Computer Output section that particularly strong corre- 
lations seem to exist between iron and copper, silver and 
copper, and magnesium and iron. 

Thus, it appears that, in general, the readings of 
different elements are not independent, and some explana- 
tion, for example, of an erroneous copper reading may come 
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from an examination of the corresponding iron reading on 
the same sample. 
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IV. ANALYSIS OF OPERATIONAL DATA 



Using the evidence provided by the Air Force data to 
support the assumption of normality in the distribution 
of the measurement error, an investigation of some opera- 
tional data was conducted. A statistical model, which makes 
use of the apparent correlations between the readings on 
different elements, was formulated. The details of the 
model's formulation are discussed in this section. 

A . NOAP DATA 

The Naval Air Rework Facility at Pensacola provided a 
magnetic data tape containing the records of operational 
analyses performed there from the beginning of July to the 
end of September in 1967. The records of some 21,000 
different analyses were included on the tape. For each 
analysis the engine model number and serial number, as well 
as the date the analysis was performed and its results in 
parts per million for each relevant element, are recorded. 

In addition, it includes the operating hours since the last 
oil change and since the last overhaul of the engine for 
each sample. Unfortunately, the action recommended by the 
laboratory after each sample was analyzed and the results 
of that action were not available with the tape. For this 
reason, there was no way of determining with certainty 
which analyses were on oil from properly operating engines 
and which were from discrepant engines. 
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Of the 113 different engine models represented on the 
tape, the Wright reciprocating engine model, R1820-82, was 
selected for investigation since it was the most frequently 
sampled, with 4,134 different analyses. 

B. PRELIMINARY RESULTS 

Since it seemed logical to expect the content of 
metallic contamination to show some increase from normal 
wear in a properly operating engine as the hours since the 
last oil change increase, some preliminary plots were made 
by the computer. Six hundred different analyses were used 
with no regard to the particular engine of the R1820-82 
type from which they came. For each of seven elements, 
relevant to the monitoring of engines of this type, the 
computer plotted the ppm content versus the operating hours 
since the last oil change. For three of the elements, 
iron, copper and aluminum, there was some indication that a 
buildup of contamination occurs. The other four plots gave 
no evidence of any significant trend. For comparison, the 
plot of chromium is included with those of aluminum, iron, 
and copper in Figures 1 to 4 , respectively. Since these 
plots were made on sample readings from a variety of en- 
gines, they did not indicate whether a particular engine 
can be expected to show the same trend. For this reason, 
the five most frequently sampled engines of the R1820-82 
type were selected and the computer was again used to plot 
the data from these engines. Five different symbols. 
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x, + , A, Q , and □ , were used on each plot to represent 
the five different engines. Thus, both the behavior of the 
readings for a given engine and differences among the 
readings of the engines could be seen. Again, evidence 
of trends was limited to the three elements iron, copper 
and aluminum. For these elements each of the five engines 
showed a roughly linear increase in the ppm content as the 
hours since the last oil change increased. For comparison, 
the plot for chromium is again included with the plots of 
aluminum, iron and copper in Figures 5 through 8. 

C. REGRESSION MODEL 

Because of the evidence of a linear increase in ppm 
content versus an increase in hours since oil change, 
provided by the computer plots, a regression model was 
suggested. With this type of model the expected content of 
a metallic element would change as the hours since the last 
oil change varies. Thus, differences between what is a 
normal amount of contamination for a properly operating 
engine just after its oil has been changed and several 
flying hours later could automatically be incorporated 
into a classification criteria. 

1 . Data Selection 

Of all the engines of the R1820-82 model type 
represented on the tape, those with eight or more different 
analyses were selected. Of these, any with missing data on 
one or more records, which brought the usable number of 
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readings to less than eight, was rejected. The remaining 
engines were then screened in an attempt to insure that all 
the data had come from properly operating engines. As was 
mentioned earlier, the tape did not provide data on either 
the recommendation made by the laboratory or the results 
of any maintenance which might have been recommended. For 
this reason, it seemed that the best way to insure that none 
of the data had come from a discrepant engine was to exclude 
all readings which exceeded the threshold limits presently 
used. These limits are given in Ref. 1. After this pro- 
cedure had been applied to the data, any engine with less 
than eight readings was eliminated from further examination. 
At this point there were 27 engines of the R1820-82 type 
remaining, each of which was represented with from eight 
to fourteen readings. This was the data that was used in 
the remainder of the investigation. 

2 . Estimation and Tests of Regression Coefficients 

It was assumed that the outcome of each spectro- 
metric analysis on a given engine is of the form 

Y = BX + e 

where Y is a seven-component vector of ppm metallic 
contents; B is a 7x2 matrix of unknown coefficients; X is 
a two-component vector with first component identically 
one and second component equal to the operating hours since 
the last oil change; and e is a multivariate normal error 
vector with mean vector zero and unknown covariance matrix. 
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Z. It will be recalled that the examination of the Air 
Force data in section III indicated that the covariance 
matrix may vary as the mean content vector changes. How- 
ever, the slopes of the regression lines mentioned in part 
C of that section are such that variations in the mean 
vectors of the extent present in the operational data do 
not result in an appreciable variation in the elements of 
£. For this reason, it will be assumed that the covariance 
matrix, Z_, is constant in the development to follow. Under 
this assumption, the matrix B associated with each engine 
was estimated using the least-squares estimation technique 
described in Appendix B. 

Since the random error vector is assumed to be 
N (0_, Z) , each observation from a particular engine is 
N (BX, Z ) . If the matrix B is partitioned so that 

B = (B 1 , B 2 ) , 

then a test of the value of each of the components of B 2 
can be made to determine whether the variability of the 
readings for a specific element is related to variations 
in the hours since oil change. The details of this test 
are included in Appendix B. For each engine and each 
component of B^ a test of the hypothesis that the component 
is equal to zero was made. Table 1 gives the number of 
times a specific component was significantly positive or 
negative at an over-all a-level of .10. For this a value, 
the expected number of times in 27 tests the results will 



32 



be significantly positive or negative, if the component is 
actually zero, is 1.35, 



TABLE 1 

NUMBER OF SIGNIFICANT REGRESSION SLOPES 
TWO-TAILED TEST; a = . 10 





Element 


Number of 
significantly 
positive 
components 


Number of 
significantly 
negative 
components 


B 2,l 


Aluminum 


8 


0 


B 2,2 


Iron 


17 


0 


B 2,3 


Chromium 


3 


1 


B 2,4 


Silver 


1 


3 


B 2,5 


Copper 


16 


1 


B 2,6 


Magnesium 


5 


0 


B 2,7 


Nickel 


5 


0 



These results indicate that the metallic content of 
aluminum, iron and copper in properly operating engines of 
the R1820-82 type tends to increase as the hours since the 
engine's oil was last changed increase. The evidence 
pointing to this conclusion is particularly strong in the 
case of iron and copper. Further, there seems to be no 
significant indication that such a relationship exists 
in general for chromium, silver, magnesium or nickel. The 
numerical results of the tests of the components of B^ f° r 
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each engine, in addition to the raw data and estimates of 
B and E_, are included in the Computer Output section. 

Any use of these results in establishing an 
operational classification procedure for identifying dis- 
crepant engines depends upon the estimation of the unknown 
matrix, B, from past analyses. For this reason, it is 
important to determine if the observations from different 
engines all come from the same over-all probability 
distribution. If this is the case, all data on engines of 
the R1820-82 model type could be used to estimate a single 
matrix, B. As a first step in this direction, the model 



Y. = B.X + e. 
— 1 —l— — 1 



was used where is a seven-component vector of readings 
on engine i? B^ = (B^, B^)^, a 7x2 matrix of coefficients 

associated with the i— engine and where the components of 
B^ associated with chromium, silver, magnesium and nickel 
are assumed to be zero; and e. is N(0, E.). The elements 

—l — —l 

of B . , not assumed to be zero, were estimated as before 

—l 

and used to estimate the 27 covariance matrices, E. . The 

—l 

unbiased estimate of E . is 



n i 

E . = Kr E (Y. . - B.X.) (Y. . - B.X.) 

-l n i -2 , =1 -1,3 -i-3 -1,3 -i-3 



. th A 

where Y. .is the j — observation of the vector Y.; B. is 
-1,3 J -i -i 

the estimate of B. ; X. is the observation of X; and 

— l — j J — 

n^ is the number of observations associated with the i— - 
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engine. A test of the hypothesis of equal covariance 
matrices was made using these estimates. The test statistic 
was extremely significant at the .10 level, and the hypo- 
thesis of equal covariance matrices was rejected. The test 
used and numerical results are included in Appendix A. 

Since the evidence indicates that the covariance 
matrices associated with readings from different engines 
are not the same, the overall conjecture of like distri- 
butions must also be rejected. 
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V. CONCLUSION 



The investigation of the Air Force data and the actual 
analysis records of a three month period from the Pensacola 
laboratory lead to three main conclusions. First, the error 
inherent in the ppm spectrometric readings is multivariate 
normally distributed with significant covariances existing 
between the readings of various pairs of elements. Further, 
there appears to be a linear increase in the content of 
aluminum, copper and iron present in properly operating 
engines of the R1820-82 type as the hours since the last 
oil change increase. Finally, there seems to be no justi- 
fication for expecting readings on samples from different 
engines of the R1820-82 type to vary in the same manner. 
Based on these results, an objective classification 
criterion can be formulated which may be of use in imporv- 
ing the present classification procedure. 

For example, all back data on a particular engine of 
the R1820-82 type which was in proper working order could 
be used to estimate the matrix B, which in turn could be 
used to estimate the covariance matrix Z_. Then, any 
observation, Y, of the spectrometric analysis of a new oil 
sample from that engine is distributed as N (BX, Z_) if the 
engine is operating properly. The estimates of B and 
could be used to construct a confidence region R (x) 
[Appendix B] . The region R^(x) would be constructed so 
that, if the engine is operating properly, the reading Y 
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will be contained in the region R^(x) with probability 1-a. 
Thus, one classification criterion would be: classify the 

engine as operating properly if Y is within R (x) and 
classify as discrepant otherwise. By making a small, say 
.01, the number of operational engines, which are mistakenly 
classified as discrepant, can be expected to be of the order 
of 1 in 100. However, the smaller the parameter a is made, 
the larger the region R (x) becomes, and thus, the more 
likely it is that a discrepant engine will be classified 
as operating properly. 

For this reason, it may be more appropriate to use 
two values of a. For example, a could be set at .10 and a' 
at .01 and two regions R (x) and R , (x) constructed. In 
this way a procedure could be used which would 1) classify 
the engine as in proper working order if Y is in R (x) ; 

OL 

2) classify as discrepant if Y is not in R , (x) ; and 3) 

" i ot 

require verification of Y or more frequent sampling if Y 

is in R , (x) but not in R (x) . 

Ot ot 

The final selection of a specific classification 
criterion and the setting of the appropriate level (s) of 
a must be done subjectively and should be based upon an 
examination of the costs involved. If the cost of classi- 
fying a discrepant engine as operational is much larger 
than the cost of grounding an operational aircraft then a 
should be made appropriately large compared to its value 
if the reverse were true. 
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APPENDIX A 



STATISTICAL TESTS 

A. TRANSFORMATION OF N (y , Z) TO N (0_, I) 

If X is a multivariate normal random variable with 
mean vector, y_, and covariance matrix, then 

Z = P(X - y) 

is multivariate normally distributed with mean vector, 0_, 
and covariance matrix, PZP 1 [Ref. 8], In addition, since 
is the symmetric matrix of a positive definite quadratic 
form, there exists an orthogonal matrix, B, such that 

BZB 1 = D 

where D is a diagonal matrix with all diagonal elements 
positive [Ref. 9]. The matrix, B, can be constructed using 
the characteristic vectors of as columns of B. Further, 
if the matrix C is defined as the diagonal matrix which has 
diagonal elements equal to the inverse of the square-root 
of the corresponding element of D, then 

CBZB'C' = CPC 1 = I 

where I_ is the identity matrix. Thus, if we define 
P = CB 



then 



Z = P ( Y - y) 

is multivariate normal with zero mean vector and identity- 
covariance matrix, and the elements of the vector Z_ are 
mutually stochastically independent standard normal random 
variables . 

B. KOLMOGOROV- SMIRNOV TEST OF GOODNESS OF FIT 

Let F(x) be defined as the cumulative distribution 
function of the random variable X which is N(0,1). In 
addition, define S n (x) to be the sample cumulative distri- 
bution function based on a set of n observations of a random 
variable assumed to be N(0,1). Then 

S n (x) = k/n 

where k is the number of observations in the sample which 
are less than or equal to x. Then, the Kolmogorov-Smirnov 
test statistic D is defined as 

D = max|F(x) - S (x) I . 
x n 1 

Observed values of D can be compared with its tabled distri- 
bution to determine the acceptability of the normal hypo- 
thesis [Ref. 10]. 

C. TEST FOR EQUALITY OF SEVERAL COVARIANCE MATRICES 
Suppose Y ^ , Y £ , , Y^ are p-component multivariate 

normal random variables with distribution denoted N(y^,Z^), 
i = l,2,...,q. In order to test the hypothesis 
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Z ^ — Z j for ci J L 1 ^ j l,2,*««,q, 

based on q samples of size I\L from the distribution of , 



l 1 f 2 f • • • 


,q, let the following quantities be defined: 


n . 
i 


= ~ 1/ i = 1 / 2 , . . . ,q ; 


n = 


q 

Z n. ; 
i=l 1 


A. 

—i 


n . 

= Z (Ki, k - - Zj.) ■ i = <3- 

k=l 



and 



A = 


q 

= Z A. . 
i-1" 1 



Then the test statistic is 



w = 


q 

= k [ Z n . (log | Z | - log | Z . | ) ] 
i=l 1 



where 



k = 


= ! _r 1 1-1] 2P 2 + 3p - 1 

i=l n i 11 6 ( P +1) (q_1) ' 

Z = A/n 



and 



/■s 

Z . 

—l 


= A i /n i i = 1,2, ... ,q. 


Asymptotic 


expansion of the distribution of this test 



statistic results in a distribution described by the 
following probability statement for W, 
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Pr (W<w) = Pr(x^<w) - c [Pr ( x^ +4 <.w) - Pr(x^<w)] - 0 (n 3 ) 
where 

*3 1 2. 2 

p(p+l) [ (P-1) (P+2) ( E -±---^)-6(q-l) (1-k) ^ ] 

i=lnf n 
i 

C = 2 

4 8 k 

f = h(q-l) (p+l)p 
2 

and Xf denotes a chi-square random variable with f degrees 
of freedom [Ref. 8] . 

This test was applied to the Air Force data as described 
in III C, with the resulting values of W = 478.5, f = 252 
and c = 5.75. The test statistic is extremely significant 
at the .01 level and the hypothesis was rejected. 

In addition, the test was used on the operational data 
as described in IV C 2 where now 

Z = A/ (n-q) 

and 

i i = A i /(n i -l) i = 1,2,. . . ,q 

instead of as defined above. The results for this test 
were W = 936.2, f = 728 and c = 14.9. Once again the test 
statistic is extremely significant at the .01 level and the 
hypothesis was rejected. 
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D. TEST FOR INDEPENDENCE AMONG A SET OF NORMAL VARIATES 
If the p-component vector Y has a multivariate normal 
distribution described by N (y_, _£) , then the test of zero 
covariance , 



E[(y i - y ± ) (y j - Pj) ] = 0, 

for all i ^ j and i,j = l,2,...,p is equivalent to a test 
of the stochastic independence of the components of Y. 
Suppose a sample, Y^ 'X.2 ' * * * '— n • n observations from the 
distribution of Y is obtained, then let V equal the deter- 
minant of the correlation matrix, R, defined by equation (2). 
The asymptotic expansion of the distribution of V results 
in the probability statement [Ref. 8] / 



Pr (-mlogV<_v) = Pr(x^<v) + -y [Pr (Xf +4 ^ v ) - Pr (x^v) ] + 0 (m 3 ) 

m 

where 

2p + 11 

m = n - , 

6 

f = hp (p-l) 

and 

n = P(P-D (2p 2 - 2p - 13) 

288 



This test was applied to the data for the Pensacola lab- 
oratory accumulated in the Air Force experiment as des- 
cribed in section III D. Table 2 gives the numerical 
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results for the ten groups of samples tested and the test 
outcomes for a level of significance of .10. 



TABLE 2 

RESULTS OF TEST FOR INDEPENDENCE 
f=21 m = 5.833 

c = 10 . 35 a = . 10 



Sample Group 


Test Statistic 


Result 


1 


26.20 


Accept 


2 


39.10 


Reject 


3 


30.81 


Reject 


4 


31.57 


Reject 


5 


38.07 


Reject 


6 


36.36 


Reject 


7 


23.78 


Accept 


8 


41.60 


Reject 


9 


23.63 


Accept 


10 


35.17 


Reject 



43 



APPENDIX B 



REGRESSION; ESTIMATION , TESTS AND PREDICTION 



A. ESTIMATION OF REGRESSION PARAMETERS 

Let Y be a seven-component random vector with distri- 
bution N (BX, E_) , where B is an unknown 7x2 matrix of 
coefficients, X is of the form (l,x) ' , and E_ is an unknown 
covariance matrix which is constant for all values of X. 

Then, the unbiased estimate of B [Ref. 8], based on a 

th 

sample of size n with the i — from N (BX^ , Z) is 

-*• — I 

B = CA 



where 



n 

C = E Y < X . 
- i=l - 1 " 1 



and 



n 

A = E X.X! 
i=l - 1 - 1 



This estimate is normally distributed with mean matrix, B, 
and covariance matrix E_ ® A ^ , where the symbol, ®, denotes 
the Kronecker product [Ref. 11]. Further, the unbiased 
estimate of E [Ref. 8] is 



E 



n-2 



n 

E 

i=l 




BX. ) (Y. 
1 — 1 



BX. ) ' 
1 



and (n-2)E_ has a Wishart distribution with parameters E^ 
and n-2. 
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B. TEST OF THE VALUE OF REGRESSION COEFFICIENTS 



Suppose the 7x2 matrix B is partitioned in the form 
B ( B , B_2 ) • 

Then for any non-null vector of constants, C, the hypothesis 
C'B = C'B* 

can be tested using an F statistic [Ref. 12] . By choosing 

t h 

C to be the vector with 1 as the i— component and all other 
components zero, and the vector B* to be the null vector, 
the hypotheses 

= 0 i = 1 , 2 , . . . , 7 

can be tested. With this selection, the standard t-test 
of the slope of a regression line [Ref. 6 and 7] can be 
used. In this way, each of the components of can be 
tested individually, each with an assigned level of 
significance, a. 

C. CONSTRUCTION OF THE REGION, R (x) 

a 

/\ A 

If the matrices B and £_ are estimated with a sample of 
size n in the manner described above, and if Y is a new 
observation from N (BX* , Z_) , the confidence region R (x) can 
be constructed so that Y will be in R a ( x ) with probability 
at least 1-a. The estimate of the mean of Y is 

A A Oi 

Y = BX 
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and has covariance matrix, 



S = T (A ® 




-1 



T' 



where 



n 

A = Z X.X! 

- i=l - 1 " 1 



and 



T = X* ® I ? . 

Thus , 

(Y - BX*) ' S" 1 (Y - BX*) 

2 

has Hotelling's T distribution [Refs. 8 and 5]. Hence, the 
set of vectors m satisfying 

(Y - m) ' S -1 (Y - m) < T 2 (a) 

comprise a 100 (l-a)% confidence region for BX * [Refs. 8 and 
5]. Since the region, R (x) , is to place bounds on Y, which 

0 i 

has covariance matrix, Z_, it can be defined by 
(Y - m) ' (S + Z) 1 (Y - m) < T 2 (a) . 



The set of all vectors, m, satisfying this constraint form 
a confidence region, R (x) , which will contain Y with 

CL 

probability at least 1 - a if Y came from N (BX* , Z_) . 
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COMPUTER CUTPUT 



A. DATA 
GRCLP A 

GRCUP A 

GRCUP A 

GRCLP A 
GRCLP A 



FRCM AIRFCRCE EXPERIMENT 



1 



AL 


FE 


CR 


AG 


CU 


MG 


N I 


3 


76 


12 


4 


62 


■a 


5 


I 


72 


8 


3 


58 


3 


1 


C 


79 


6 


2 


63 


5 


1 


1 


79 


10 


4 


61 


3 


5 


1 


12 


4 


2 


61 


2 


1 


G 


75 


9 


4 


6 1 


a 


1 


3 


73 


1 1 


4 


59 


3 


2 


1 


IS 


3 


2 


66 


3 


l 


0 


74 


2 


3 


62 


3 


0 


AL 


FE 


CR 


AG 


CU 


MG 


NI 


c 


127 


7 


1 


12 


10 


2 


C 


142 


1C 


2 


13 


1 1 


3 


1 


144 


15 


3 


14 


1 1 


3 


0 


147 


9 


2 


14 


1 1 


3 


c 


143 


9 


3 


14 


11 


3 


c 


154 


14 




15 


12 


4 


r 


145 


17 


4 


15 


1 1 


4 


C 


153 


14 


3 


15 


13 


4 


2 


77 


1 1 


5 


6 


14 


1 


C 


84 


8 


4 


8 


16 


2 


AL 


FE 


CP 


AG 


CU 


MC- 


N I 


C 


15 1 


14 


3 


16 


12 


3 


I 


151 


13 


2 


14 


1 1 


4 


1 


84 


9 


4 


8 


16 


2 


c 


84 


12 


4 


6 


15 


l 


c 


ec 


1 C 


4 


c 


14 


2 


c 


82 


7 


4 


7 


15 


1 


c 


83 


11 


5 


8 


16 


1 


c 


88 


6 


4 


7 


17 


0 


I 


77 


7 


4 


6 


14 


1 


1 


82 


8 


3 


8 


15 


1 


AL 


FE 


CR 


AG 


CU 


MG 


N I 


26 


27 


9 


4 


35 


1C 


13 


23 


27 


8 


4 


36 


1 1 


13 


24 


27 


15 


4 


37 


1C 


14 


20 


28 


11 


4 


34 


1 1 


13 


24 


28 


15 


5 


34 


1C 


14 


23 


2 9 


12 


5 


34 


1 1 


14 


22 


28 


10 


5 


34 


1C 


13 


23 


28 


14 


5 


35 


1C 


14 


AL 


FE 


CR 


AG 


CU 


MG 


NI 


2 


137 


9 


4 


ICO 


6 


1 


5 


132 


17 


6 


107 


c 


1 


6 


127 


17 


6 


1C3 


5 


1 


3 


133 


2 


4 


116 


G 


1 


4 


141 


4 


3 


111 


6 


0 


■a 


138 


6 


4 


117 


6 


1 


2 


135 


12 


5 


1 C 3 


c 


1 


4 


145 


18 


6 


1 14 


6 


2 


3 


137 


1 1 


5 


1C9 


5 


1 


3 


133 


1C 


5 


1 C 4 


6 


1 
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GRCUP # 



GPCUF * 



GRCUP * 



GRCUP # 



GRCUP * 



AL 


F F 


CP 


AG 


CU 




MI 


r 


112 


11 


7 


13 


19 


4 


0 


1 ") 5 


c 

-r' 


2 


11 


15 


2 


1 


1C. 7 


9 


3 


13 


1 9 


2 


2 


Q Q 


10 


4 


13 


17 


3 


C 


1C 9 


8 


2 


12 


18 


2 


C 


1 ? 4 


11 


3 


15 


1 c 

^ • 


3 


r*. 


lie 


1C 


3 


1* 


20 


4 


r 


115 


11 


3 


15 


21 


3 


i 


1 10 


13 


3 


1 4 


?c 




i 


125 


11 


3 


16 


22 


3 



A L 


F E 


CR 


AG 


CU 


W G 


MI 


3 


95 


4 


4 


85 


27 


1 




95 


8 


4 


70 


7 K 


1 


1 


87 


10 


2 


72 


25 


2 


3 


102 


11 


4 


86 


27 


1 


1 


1 94 


10 


4 


90 


2 8 




2 


1C 


10 


5 


70 


27 


2 




105 


11 


R 


8 3 


28 


5 


£ 


l n 6 


11 


4 


37 


23 


1 


i 


1C 4 


8 


4 


81 


27 


l 



A L 


F E 


CR 


AG 


CU 


M G 


MI 


2 


132 


6 


3 


5 


11 


1 


r 


122 


8 


5 


c 

> 


11 


3 


r 


127 


4 


3 


8 


11 


3 


r 


133 


11 


4 


-* 

r 


12 


4 


1 


14 D 


4 


7 


6 


12 


2 


1 


138 


1 


2 


6 


13 


1 


f 


127 


2 


3 


6 


11 


2 


r 


13C 


2 


7 


8 


12 


2 


r 


137 


7 


4 


7 


1 7 


3 


r 


123 


9 


4 


8 


I 1 


4 


C 

AL 


FE 


CR 


AG 


CU 


MG 


MI 


2 


17 


12 


7 


9 


12 


1 


1 


1 3 


1 1 


2 


10 


12 


0 


7 


1 7 


11 


2 


n 


1 7 


1 


c 


1 7 


3 


1 


1C 


i? 


0 


1 


17 


13 


7 


10 


12 


0 


1 


19 


ir 


2 


11 


13 


1 


I 


1 3 


6 


1 


11 


12 


1 


1 


1 8 


9 


2 


1 1 


12 


1 


AL 


F E 


CP 


AG 


CU 


M G 


MI 


r 


187 


8 


3 


7 


10 


4 


4 


183 


10 


3 


8 


10 


7 


1 


197 


12 


3 


o 


9 


4 


n 


19 0 


8 


3 


7 


1 A 


4 


c 


j 


0 


o 


o 


n 


C 


r 


188 


8 


2 


Q 


10 


4 


i 


212 


7 


2 


10 


10 


3 


r 


206 


4 


1 


10 


11 


4 


C 


217 


8 


2 


10 


11 


5 
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B 



SAMPLE MEANS ANC STANCARC CEVIATIGNS 



GROUP 


ft 


l: 


MEAN 

S.D. 


1 


AL 

1. 11 
. 167 


GROUP 


ft 


2: 


MEAN 

S.C. 


0 


AL 

0.20 

.422 


GROUP 


ft 


3 : 


MEAN 

S.D. 


C 


AL 

C .50 

. 7C 7 


GRCLP 


ft 


4 : 


MEAN 


AL 

23. 12 


S.D. 


1 


.727 


GRCLP 


ft 


a • 

v * 


MEAN 

S.D. 


1 


AL 

3.50 

.265 


GROUP 


ft 


6: 


MEAN 

S.D. 


C 


AL 

0.50 

.707 


GRCLP 


ft 


7: 


MEAN 

S.C. 


l 


AL 

3.11 

.537 


GROUP 


ft 


8: 


MEAN 

S.C. 


0 


AL 

C .40 
.699 


GROUP 


ft 


9: 


MEAN 

S.C. 


0 


AL 

1.25 

. P 86 


GROUP 


ft 


10: 


MEAN 

S.C. 


o 

J 


AL 

C. 33 
.500 



EE 


CR 


75.56 


7.22 


2.835 


3.632 



FF 


CR 


145.70 


12.20 


7.8^8 


3.225 



EL 


CR 


8 2.10 


8.90 


3.384 


2.C25 



FE 


CR 


27.75 


11.75 


C.7C7 


2.712 



FE 


CR 


136.20 


1C .60 


5.119 


5.582 



FE 


CR 


1C9. 70 


9 . 5 C 


7.530 


2.183 



FE 


CR 


ICO. 22 


9.22 


6.039 


2.279 



FE 


CR 


131.40 


5. 40 


5.702 


3 .340 



FF 


CR 


17.62 


10. CO 


C .744 


2 .268 



FE 


CP 


158.56 


8. CO 


11.294 


2.062 



AG CL 

3.11 61.44 

C.528 2.299 



AG CU 

2.60 14.20 

0.843 1.135 



AG CL 

4.10 6.50 

C . 5 6 8 l .101 



AG CU 

4.50 34.68 

0.535 1.126 



AC CL 

4 . 6 C 108.40 
1.033 5.969 



AC- CU 

2.80 13.70 

0.632 1.567 



AC- CU 

4.CC 82.78 
C.866 5.653 



AG CU 

3.20 6 . 6 C 

1.033 1.174 



AG CU 

1.67 10.25 

C.641 0.7C7 



AG CU 

2.33 8.78 

C.7C7 1.2C2 



MG M 

2.89 1.85 

0.333 1.833 



MG M 

11.30 3.30 

0.823 0.675 



MG M 

15.20 1.20 

1.033 0.632 



MG M 

1C. 37 13.50 

0.518 C.535 



MG N I 

5.50 l.CC 

0.527 C.471 



MG M 

15.40 2.90 

1.430 C.738 



MG N I 

26.89 1.67 

1.168 1.323 



MG N I 

11.60 2.50 

C.655 1.C8C 



MG N I 

12.25 C . 6 3 

0.463 0.518 



MG N I 

10. 11 4.33 

0.601 1.118 
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SAMPLE CORRELATION MATRICES 



GRCtP ft l: 
AL 

AL l.COO 
FE -C.21C 
CR 0.6 A? 
AG C.449 
CU -0.254 
MG C .036 
N I 0.591 



r E CR 

-C.210 0.642 

l.COC -C.C50 
-0.050 1.000 

-C.169 0.808 

C.746 -0.5C7 

0.338 0.333 

C.35C C.736 



AG CU 

0.449 -0.254 

-C.169 0.746 

0. 8C8 -C.5C7 

1. CCC -0.495 
-0.495 i.OCC 

0 .449 0 .073 

0.596 -0.076 



MG M 

0.C36 C.591 
C .338 C.350 
0.333 C.736 
0.449 C .596 
0.073 -C.C76 
i.OOC C.182 
0.182 1.000 



GRCUP 


ft 2: 
AL 


FE 


CR 


AL 


1 .000 


C . 12 1 


0.294 


FE 


0.121 


1 .000 


0 .631 


CR 


C.294 


0.631 


1 .GCO 


AG 


-C .062 


0.568 


C.8C9 


CU 


-C.093 


C.831 


0.747 


MG 


-C . 192 


C. 807 


C.519 


N I 


0 .156 


0.816 


0.735 



AG 


CU 


MG 


Nl 


-C.C62 


-0.093 


-C. 192 


C.156 


0.568 


0 .831 


0 .807 


C.816 


C .809 


0.747 


0.5 19 


0 .735 


l .000 


0.789 


C.512 


C.625 


0.789 


1 .000 


0.761 


C .636 


C.512 


0.761 


1 .000 


C .620 


C.625 


0 .638 


C .620 


1.000 



GRCUP 


ft 3: 
AL 


FE 


CR 


AG 


CU 


MG 


NI 


AL 


I.OCC 


-C.627 


0.116 


0.138 


-0.071 


-0 .456 


0 .000 


FE 


-0 .627 


1.000 


-C.225 


-C.237 


C.510 


0.916 


-0.218 


CP 


0.116 


-0.225 


l.OCO 


C.493 


-0 .254 


-C .308 


0.278 


AG 


0.138 


-C.237 


0.493 


1 .000 


-0 .160 


-0 .C38 


-0 .062 


CU 


-C .071 


C. 510 


-0.254 


-C .16C 


l.CCC 


C.704 


C.C32 


MG 


-0.456 


C.916 


-0 .308 


-0 .038 


0 .7 C4 


1 .000 


-0.238 


NI 


c.coo 


-0.218 


0.278 


-0.062 


0.0 32 


-0 .238 


1 .000 



GROUP 


ft 4 : 
AL 


FF 


CR 


AG 


CU 


MG 


NI 


AL 


1.000 


-C.439 


o.ceo 


-C .077 


0.377 


-C.54C 


C .232 


FE 


-0.439 


1 .000 


0.261 


0 .756 


-0.763 


0.293 


0.378 


CR 


C .069 


C.261 


l.CCO 


C.354 


0.082 


-0.433 


C.887 


AG 


-0 .077 


C.756 


0.394 


I.OOC 


-0.593 


-0.258 


C .50C 


CU 


C.377 


-C.763 


0.082 


-0 .593 


1 .000 


-0.153 


0.119 


MG 


-0.54C 


0.293 


-0.433 


-C.258 


-0.153 


l.CCC 


-0.258 


NI 


C .232 


0.378 


0.887 


C .50C 


0.119 


-0.258 


l.COO 



GRCUP ft 5: 





AL 


FE 


CR 


AG 


CU 


MG 


NI 


AL 


l.COC 


-C.444 


0.518 


0.509 


0.044 


-0 .249 


C .000 


FE 


-C .444 


1.000 


-0. 110 


-0.286 


0.369 


C.536 


C.184 


CR 


0.518 


-C.liC 


1.000 


C .929 


-0.385 


-0.227 


0.591 


AG 


C.5C9 


-0.286 


0.929 


l.OCO 


-0.256 


-0 .408 


0.685 


CU 


0 .044 


C .369 


-0.385 


-C .256 


l.OCO 


0.142 


C.115 


MG 


-C . 249 


C.536 


-0.227 


-0 .408 


0.142 


1 .OOC 


C .000 


NI 


C.OCC 


C. 184 


0.591 


0.685 


0.119 


0.000 


1 .000 
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GROUP « 6 



AL 
FE 
CF 
AG 
CU 
MG 
N I 



GROUP 



AL 
FE 
CP 
AG 
CU 
MG 
N I 



GROUP 



AL 
FE 
CR 
AG 
CU 
MG 
N I 



GROUP 



AL 
PE 
CR 
AG 
CU 
MG 
N I 



GROUP 



AL 

FE 

CR 

AG 

CU 

mg 

N I 



AL 

1.000 
- J. 198 

3.252 

0.74* 

3.0 50 
-0. 7 20 
- 3 • 106 



H 7: 
AL 

l.oro 

0.469 
3.206 
1. 376 
3.593 
j. 426 
-7.164 



H 3: 

AL 

1.000 

0.457 

-3.219 

-7.431 

-0.596 

•0.136 

-3.736 



H 9: 
AL 

l.cco 
-7.271 
7.426 
3.566 
-7. 342 
7.522 
3. 5 h5 



H 10 : 
AL 

l.OCO 
7.C3C 
0 • 5 C 6 
7. 354 
3.139 
-7. 555 
0.224 



Ff 

-C .198 
l.C^'C 
C. 322 
-C.2CC 
C. * 75 
0. 89 f 
C. 2 14 



FE 

C • 4 6 9 

1. 30C 
C • 2 0 5 
0.741 
C.764 
C . <*(. c 
t.C 89 



FE 

C.4f7 
1 • C < C 
-C . l c 6 

-c. 543 

-C.U-4C 

C.77C 

-C.343 



FP 

-0. 271 
1 . -• C 0 

-0.339 
-C. 112 
c. ^47 

C. 311 
C.325 



PF 

f . 03C 
1 . - or 
-C. 376 
-C.621 
C. 83C 
C.5C6 
-C.264 



CP 

0.252 
0. 322 
1.00 0 
0.467 
3.737 
0.370 
0.614 



CR 

0. 2C6 
0.20 5 
1.000 
0. 063 
- 0 . 122 
0.246 
0. 359 



CR 

-0.219 
-C. 196 
l.CCC 
0. 812 
0.045 
-0. 3C4 
0. 770 



CP 

7.426 
-0. 339 
1.000 
0. 7 86 
- 3. 624 
7. 136 
-0.243 



CR 

0.606 
-0. 376 
l.COO 
0. 772 
-0. 303 
-0. 7C6 
0. 380 



AG 

0.745 

- 0.200 

0.467 

1.000 

0.493 

-C.C24 

0.191 



AG 

0.3 76 
G. 7 41 
0.063 
1.0 cO 
0.460 
(-.613 
0.327 



AG 

-0.431 
-0.543 
912 
1.070 
- 0.110 
-0.492 
0. 697 



AG 

C. 566 
-C. 112 
0.786 
1.000 
-7. 552 
0 * 120 
0.269 



AG 

0. 354 
-C.621 

0. 772 

1. DOC 

-0.784 

- 0.686 

0.316 



CU 

0.^5 •'< 
575 
0. 737 
0.493 
1 . 0 0 0 
7.7C4 
7.5 48 



CIJ 

‘".593 
7. 764 
-' 1.122 
0.460 

1.0 no 

9.773 

-7.228 



CU 

-0 .596 
-0. 04 n 
7.045 
- >.110 
l.OCC 
0.0 54 
7.526 



CU 

342 

7.747 

-7.624 

-9.552 

1.000 

0.218 

0.293 



CU 

7.139 
0. 8 30 
-0.303 
-0.734 
l.OCO 
0.38* 
-0.217 



MG 

-C .220 
0.890 
0.37C 
-0.024 
0.704 
1.000 
0.253 



MG 

C • 4 2 6 
0.90© 
0.246 
0.618 
0.77 3 
1.000 
0.216 



MG 

C. 136 
C-.77C 
— 0 . 30 4 
-0.492 
0.054 
l.OCO 
-0.294 



MG 

0.522 

0.311 

0.136 

0.12C 

0.218 

1.000 

0.447 



MG 

-0.555 
7.506 
-0.7C6 
-Cl 686 
0.38 5 

1.000 

0.124 



MI 

- 0 . 1 n 6 
0.214 
0.614 
0.191 
0.548 
0.253 
1.000 



MI 

-0.164 

7.359 
0. 327 
-0.228 
0.216 

1.0 or 



NI 

-0.736 
-0. 343 
0.770 
0.697 
0.526 

-0.294 

1.000 



MI 

0. 545 
0.325 
-0.243 
0. 26° 
0. 2Q3 
0.447 
1.0 70 



NI 

0.224 

-0.264 

J.38C 

0.316 

-0.217 

0.124 

1.000 
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D. RESULTS 



ALLRI NUM 

INTERCEPT 
SLOPE 
CORPELAT ICN 



I RCN 

INTERCEPT 
SLOPE 
CORPELAT ICN 



CFRCRILM 

INTERCEPT 

SLOPE 

CORRELATION 



SILVER 

INTERCEPT 

SLOPE 

CORRELATION 



COPPER 

INTERCEPT 

SLOPE 

CORRELATION 



RAGNES. 

INTERCEPT 
SLOPE 
CCRRELAT ICN 



NICKEL 

INTERCEPT 
SLOPE 
CGRRELAT ICN 



RESULTS OF 
REGRESSION 



OF REGRESSION CF VARIANCE CN -H E AN SQUARED 



0 . 8768 E OC 
0 .4C16E-C2 
C . 7C6?E CC 



RESIDUAL VARIANCE 
VARIANCE OF INTERCEPT; 
VARIANCE OF SLCPE 



C.4118E Cl 
0.5782E-C1 
C. 1824E-C5 



•0.4542E Cl 
0 .3C8CE-C2 
0.9296E OC 



RESIDUAL VARIANCE 
VARIANCE OF INTERCEPT 
VARIANCE CF SLCPE 



0.1829E 04 
C.5537E 02 
C.83C5E-C7 



C.7276E 01 
0 .2642E-01 
0.1 1 9 3 E OC 



RESIDUAL VARIANCE 
VARIANCE OF INTERCEPT 
VARIANCE CF SLCPE 



C.5985E C3 
C.5717E C2 
C • 8870E-C3 



0.5494E OC 
0 .7357E-C2 
0 . 16 5 E CC 



RESICUAL VARIANCE 
VARIANCE OF INTERCEPT 
VARIANCE OF SLCPE 



0.7648E 00 
0 .4397E-01 
C.58C0E-04 



0.3822E CC 
0 .3 238E-C2 
0 • 9 4 94 E OC 



RESIDUAL VARIANCE 
VARIANCE OF INTERCEPT 
VARIANCE CF SLCPE: 



C.1631E 03 
C.2883E Cl 
0 . 1054E-C6 



: 0.2444E CC 
: 0.2224E-02 
: 0.7544E CC 



RESICUAL VARIANCE: C.15C3E Cl 
VARIANCE CF INTERCEPT: 0.3734E-C1 
VARIANCE CF SLCPE: C.2511E-C6 



: 0.1C66E Cl 
:-0 .4177E-C2 
:-0.237CE OC 



RESICUAL VARIANCE: C.8284E Cl 
VARIANCE OF INTERCEPT: C.1239E 00 
VARIANCE OF SLOPE: C.3126E-C4 



A T TEST CF THE HYPOTHESIS THAT THE SLCPE OF THE 
LINE IS ECUAL TO ZERO 



ALUMNUN 


0.2973542E 


01 


REJECT 


IRON 


0. IC68822E 


C2 


REJECT 


CHRCMIUN 


0 . 886958 1 E 


00 


ACCEPT 


SILVER 


C.9660166E 


00 


ACCEPT 


CCP PER 


0.9976842E 


01 


REJECT 


M AGNES . 


0 .4439077E 


01 


REJECT 


NICKEL 


-0.7471411E 


00 


ACCEPT 
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E. OPERATIONAL CATA, ESTIMATED PARAMETERS AND TESTS OF 
REGRESSION SLOPES FOR ENGINES OF THE R182C-82 MCCEL 

RESULTS OF LINEAR REGRESSION USING THE MODEL, 

Y = BX + E 

WHERE: 

Y IS A VECTOR OF WEAR METAL CONTAMINATION IN OIL 
SAMPLES 

Y I = PPM OF ALUMINUM 

V 2 = PPM CF IRON 

Y2 = PPM CF CHROMIUM 

>4 = PPM CF SILVER 

Y 5 = PPM CF COPPER 

Y 6 = PPM CF MAGNESIUM 
Y7 = PPM CF NICKEL 

B IS A 7X2 MATRIX CF COEFFICIENTS 

X IS A VECTOR CF VARIABLES 
XI = 1 

X 2 = HOURS SINCE OIL CHANGE 

ANC fc IS A 7 COMPONENT MULTIVARIATE NORMAL RANDOM 
ERROR VECTOR WITH MEAN VECTOR ZERO ANC UNKNOWN 
COVARIANCE MATRIX SIGMA 
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ENGINE NUMBER: 316726 



DATA : 



Y 1 


Y 2 


Y 3 


Y4 


Y5 


Y 6 


Y7 


XI 


6 . 0 


2 5.0' 


6 . U 


2. J 


12.' 


3.0 


1.0 


11.0 


E. r 


2 8.0 


8. ) 


2.0 


11 . 0 


1.0 


2.0 


57, 7 


1C. 


4 1.0 


3. 


3. 0 


2 3.0 


2 • C 


2.0 


256.0 




2 5 . 0 


l.o 


2.0 


12.0 


1 .0 


1.0 


21.9 


9.o 


31.o 


6.0 


2.o 


10.0 


2.0 


l.o 


ft A, 


1C. 


3 ') . 0 


L.o 


l.o 


I 6.0 


2.0 


o.O 


89.0 


6. •. 


26. 0 


3 . >0 


2. 0 


14.o 


1.0 


£ /> 
. 


64.0 


1 1. r 


3 ) . 9 


1 . 0 


1.0 


I 6.0 


1.0 


1 • ( 


133. A 


7.; 


27.0 


3.0 


r ■> 

v. • , 


12." 


l.o 


2.0 


l n r . 0 . 


11.0 


3 8.0 


7 . v 


1. T 


19.0 


2.0 


2 . C 


171. " 



NUMBER OF DATA PCIMS = 10 
RESULTS: 

f ST I MATE C f R 



81 


82 


7 . 1 4 7 1 . 6 


.0 15 62 


23.26 556 


0 . 0 6 8 7 5 


4. '35f° 


-0. 00,2 4 1 


1. 5C 81. 4 


O.OO 093 


5. 71871 


0. 34810 


1 . 64-72? 


6. 'MO 5 3 


1 . 9340 ? 


0.00066 



ESTIMATE CF CCVARIANCE MATRIX SIGMA 



2.54-2 1.314 

1.314 4.626 

- 1 . L 9 P 3.393 

-C . 2 3 7 J • 7 Q 2 

C . 4 r c 0.883 

-C.f-75 1.955 

-C.937 -5.65: 



-1.098 -0. 237 

3.390 9,792 

8.665 0.538 

0.538 0.795 

-2.176 0.402 

C .658 0.172 

-0.615 -j.904 



J.409 -0.075 

',883 o . 9 R 5 
- 2.176 r '.65 8 

1.482 0.172 

4.619 0.471 

r .471 ^.543 

'•. 80 8 - 0.25 2 



-u. 937 
- 0 . 650 

-0.61 



T TEST CF THE HYPOTHESIS THAT A 
COMPONENT CF THF R2 V C CTCP FCUALS 
ZE°C AT A .10 SIGNIFICANCE LEVEL 



TL ST STATISTIC RESULT 

2.15173 REJECT 

7.(1834 REJECT 

-C. 17978 ACCEPT 

C. 22784 ACCEPT 

4 . c 1 42 3 REJECT 

0.15746 ACCEPT 

f. 11025 ACCERT 



54 



-4INJ OD.tsiT\ 



ENGINE NIJMRER: 515355 



C A TA : 



Y 1 


Y2 


Y 3 


Y4 


Y 5 


Y6 


Y 7 


XI 


K . ' 


3 9 . 0 


5.0 


c.o 


17.0 


2. : 


4.C 


126.0 


9.0 


?R. 0 


2.0 


1.3 


11.0 


2 . C 


3.0 


117.0 


12. * 


31.0 


3.3 


0.3 


14." 


2.0 


2. C 


84. 0 


6.0 


2 3.0 


0.0 


2.0 


10.0 


o . C 


o.O 


66.0 


5 • • 


13.0 


5.n 


0.0 


7 . 0 


o.o 


l.C 


24.0 


4. ‘ 


2 J. r 


3.0 


l. T 


7.0 


1.0, 


2 . C 


10. o 


8. 


"3 T > I 


1 • o 


0. ) 


15.° 


2.0 


1.0 


79.0 


C 

- • 


26. > 


6 . n 


2. ' 


12.0 


1 . 0 


l.C 


46.o 


1C .r 


2 8.0 


5.0' 


3. 0 


ll.o 


2. 0 


4.0 


15.0- 



MJMPEF PE GAT 4 POINTS = 5 



PE SLIT 5: 

estimate cf e 



n 1 

5. 1 936* 
1 C .H3^9 C 
6. 56766 
1.73114 
7.73 G 36 
C.t P'42 
1.4 3Af 2 
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U . C 3 P. 9 8 

3.11819 
-0.03657 
-0. Cl 152 

’)• C6C 15 
0 • 0 1 n 2 3 
0. OCR 91 



ESTIMATE CF CCVARIANCE MATRIX SIGMA 



787 


6.406 


2.578 


0. 336 


2.786 


1.337 


2. 


013 


4f 6 


20.915 


4.30 5 


0.46 6 


1 r > . 3 G 7 


2.700 


3. 


444 


578 


4. 305 


9.461 


- 3,004 


1.62C 


0.526 


2. 


66 5 


336 


0,466 


-0.004 


1.160 


1.115 


1.095 


r- 

V * 


35 r 


786 


10. 397 


1 . 6 2 C 


0.115 


5.864 


1.090 


3 . 


774 


337 


2.70* 


0.526 


3.095 


1.C9C 


0.64 6 


0. 


816 


^13 


3.444 


2.665 


0. 350 


0.774 


9.816 


2. 


125 



T TEST CF 


THE HYPOTHESIS THAT A 


COMPONENT 


CF THE 62 VECTCR 


ECUALS 


ZE p C AT A 


.K SIGNIFICANCE 


LEVEL 


T STATISTIC 


RESULT 




1.92706 


REJECT 




3.07322 


REJECT 




-1.^ 2541 


ACCEPT 




-1.272*1 


ACCEPT 




2.95414 


REJECT 




1. 51371 


ACCEPT 




C. 72672 


ACCEPT 
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ENGINE NUMBER: 51667a 



DATA : 



VI 


Y 2 


Y 3 


Y4 


Y 5 


Y 6 


Y7 


XI 


5. S' 


25.9 


5.0 


1.0 


13.0 


3.0 


5 • C 


60.0 


P.' 


2 5.0 


7.0 


1.0 


12.0 


1 .C 


5 . n 


89.0 


P . 


2d. 0 


3.0 


1.0 


1 l.r 


1.0 


1 . o 


75.0 


7.. 


27. 7 


2.0 


0, 0 


14.0 


2 . ^ 


5.0 


5 8. 0 


1 C . ■; 


35. 0 


5.0 


0. 9 


14.0 


3. r 


3.0 


25.0 


P." 


31.C 


l. r 


0.0 


12. r 


3.U 


4.0 


3 1 . 0 


6. • 


3 •) . 0 


8 • G 


1 . 0 


12.0 


4.0 


7.0 


117.0 


6, 


26. C 


5.0 


2.0 


in.?'! 


2.0 


3.0 


51. C 


5 »<" 


2 3.0 


9.0 


1.0 


a. : 


1.0 


o.o 


18.0 



NUMBER OF TATA PCINTS = 9 

R E SUt T 5 : 

ESTIMATE CF R 



01 


P 2 


6. 73207 


0.312 0 5 


26.238? 0 


0. CO 5 37 


3.9< 117 


•o,0 179 1 


0.4 c 374 


0.00633 


1 1.0 53 23 


0 . C 1 2 2 6 


1.8587 3 


0.00615 


2.4 8 8 7 r 


0.71425 



ESTIMATE CF COVARIANCE MATRIX SIGMA 



8 64 


6. 322 


-1.535 


-3.961 


2.956 


0.Q3O 


1. 


468 


322 


27,427 


-2.357 


-1.023 


4.368 


3. 946 


0 . 


955 


535 


-2. 397 


7.515 


0. 865 


-3.111 


-0.413 


-2. 


296 


961 


-1.023 


3.869 


0.46 2 


-1.91C 


-3.267 


-0 . 


4 38 


956 


4.363 


-3.111 


-1.010 


4.048 


1.119 


2. 


893 


CD,’; 


3. 946 


-0.41 3 


-0.267 


1.119 


1.321 


0. 


8 30 


468 


0. 959 


-2.296 


-0.438 


2.893 


0.30 3 


4. 


0 5 C 



T TEST OF 


THE HYPCTHES IS THAT A 


component 


CF THE 02 VECTCP 


ECU AL S 


ZFRC AT A 


.10 SIGNIFICANCE 


LFVEL 


t statistic 


RESULT 




0.55162 


ACCEPT 




C • r 5 7 2 0 


ACCENT 




0.57286 


ACCEPT 




f .83792 


ACCEPT 




0 . 5 4 8 1 C 


ACCEPT 




5.43131 


ACCEPT 




0 .63835 


ACCEPT 
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ENGINE NUMBER: 515058 



DATA : 



Y 1 


Y2 


Y3 


Y4 


1C.. 


2 3.3 


7.0 


1. 1 


7.: 


25.0 


2.0 


0 . 0 


4.'' 


20. 0 


6.0 


0.0 




2 3.0 


•: . o 


o.O 


6. 0 


24.0 


5.'' 


l.o 


7.-. 


20.0 


7.0 


1.0 


6. r 


15.0 


3.0 


1. 0 


c * 

^ • 


26.0 


e.o 


0. -t 


5.o 


18.0 


3.0 


1.0 


NUMRER PF 


DATA 


PCINTS = 


q 


PESULTS : 









ESTIMATE CF P 



ni 


P 2 


P. c 3^6-" 


-0 • C 1 9 84 


2t ,.2^82° 


) • C 32 2 0 


7.1^537 


-O.C3448 


C. opr H 


-O.OOFf s 


6. 16235 


0 . 0 1 8 5 5 


-f . 290^7 


0.:1422 


1.21374 


0.C03O7 



Y 5 


Y 6 


Y7 


XI 


8.0 


0.0 


2.0 


2 6."? 


lO.o 


2.0 


3. n 


1 3 c . 0 , 


0.0 


1.0 


2.0 


8 3.0 


7.0 


2.6 


1.0 


73.o 


9.0 


1.0 


l.C 


107.^ 


9.0 


0 . o 


0 . 0 


sr. n 


6.0 


1.0 


3.C 


44.0 


7.0 


0.0 


l.C 


5 6.0 


4.0 


o.O 




8 8.0 



ESTIMATE CF CCVARI A N C S MATRIX SIGMA 



1 7 ♦ 5 r 9 


18.187 


- 7.369 


-0.609 


’0 .986 


2 , 7 7 3 


-c. 


173 


18. 1P2 


3 n .052 


-3.0C 2 


-1.663 


2. 588 


1.725 


-0. 


941 


-7. 36 c 


-3. "02 


6,738 


0.059 


2.0 24 


-1.333 


-r. 


60 5 


- 0 . 60c 


-1.663 


0 . 059 


0 . 275 


-0.115 


-0.163 


- 0 . 


152 


C. 8 P 6 


.3.588 


2.024 


-0.115 


3.28? 


237 


c. 


607 


2,073 


1.725 


-1.333 


-0.163 


•0.237 


0.52 5 


0 . 


40 s 


-C.173 


-0.941 


-0.605 


-0.152 


0 . 6 C 7 


0.4Q8 


1. 


448 



T TEST OF THE HYPOTHESIS THAT A 
CCVPCNENT CF THE B2 VECTOR ECUALS 
ZFRC AT A .1C SIGNIFICANCE LEVEL 



TEST STATISTIC 


RESULT 


-C. 45716 


ACCEPT 


0.56631 


ACCEPT 


-1.2 6071 


ACCEPT 


-1.03937 


ACCEPT 


C. 9«61 5 


ACCEPT 


1.89257 


ACCEPT 


0. 24613 


ACCEPT 
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ENGINE NUMBER: 509090 



rATA : 



Y 1 


Y2 


Y 3 


Y4 


21.0 


29. 0 


3,0 


0.0 


16. j ' 


2 2.0 


2.C 


1.0 


18.” 


4 ■> . o 


0 . 6 


3. 0 


17.0 


2 9.0 


0.0 


j.O 


9. o 


2 9. 0 


5 . C 


0. 0 


P.D 


21.0 


7.0 


1.0 


1 4 . ( 


3 7.0 


6.0 


Q.O 


e. •; 


26. C 


0 . C 


0. 9 


15. <• 


4 4.0 


5 • C 


1.0 


1 2 • C 


3 3.0 


3, r 


0.0 


16. 


3 3.9 


4 • C 


0. 3 


NUMPEF LF 


DATA 


PCI NTS = 


1 1 



RESULTS : 




F S 


TIMATE CF P 


PI 


92 


9.01636 


0.03755 


18.23833 


6.10673 


3. P4371 


-0.0049a 


r. 73897 


-0.00351 


4 , 1 1 8 8 z 


C. 04157 


C • 3 7 8 3 ° 


^ 5 


1 . 36356 


-0.0 0160 



Y 5 


Y 6 


Y 7 


XI 


8.0 


3.0 


l.C 


111.0 


7 .:- 


2.6 


1 . c 


90.0 


12.0 


1.0 


1.0 


1 6 1 . C 


11.6 


■l.C 


l.o 


183. 0 


7.0 


r.c 


l.C 


9C.C 


7. : 


f . 0 


1.0 


30.0 


12.0 


l.C 


• 0 


198.0 


6.6 


o.c 


l.o 


60,0 


11.6 


2.0 


4.0 


156. ? 


12.0 


l.C 


6.6 


193.0 


13. ” 


2.0 


l.C 


188.0 



C . 783 
3. 044 
C. 282 
0,281 
0, 187 
0 . 351 
1.2 0 2 



T TEST CF T H F HYPOTHESIS THAT A 
COMPONENT G F THE C2 VECTOR FC'JALS 
ZFRC AT A .10 SIGNIFICANCE LEVEL 



FSTIMATE CF COVARIANCE MATRIX SIGMA 



14, 99 1 3.956 

0.95 i 20. 278 
-3.165 -0.155 

0.176 0.297 

—C , 366 1.499 

3.122 -0.153 

C.788 3.944 



-3.165 0.174 

-0.155 0.297 

6.752 0.427 

0.427 0.195 

C.771 0.130 

-C .035 0.133 

0 . 2 8 2 0 . 281 



-0.366 3.122 

1.49 9 - r . 153 

0.771 -0.035 

0.130 6,133 

0.939 -0.114 

-0.114 0.929 

0 . 187 0 . 3 R 1 



TEST STATISTIC 


RESULT 


1.89834 


ACCEPT 


4.42040 


REJECT 


-0.35791 


ACCEPT 


-1.* 9459 


ACCEPT 


8.0)668 


REJECT 


1.17122 


ACCEPT 


-0.27252 


ACCEPT 
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ENGINE NUMBER: 524739 



DATA: 



Y 1 


Y2 


Y3 


Y4 


Y5 


Y 6 


Y7 


XI 


1 6 . C 


22.0 


8.C 


2.0 


4.0 


0.0 


2.0 


2P.0 


13..) 


37.0 


O.C 


0. o 


13.0 


2.0 


3.0 


181. C 


1 6. ° 


35.0 


1.0 


2.0 


11.0 


5.0 


6.0 


32.0 


6.o 


29.0 


O.C 


3.0 


9.0 


2.0 


3.0 


158.0 


5.C 


27.0 


1.0 


0. 0 


7.0 


O.C 


l.C 


15.0 


16.R 


39.0 


1.0 


2.0 


11.0 


5.0 


3.0 


6 2. C 


1 6. r 


*5.0 


2.0 


4. ^ 


14.0 


4.0 


1.0 


169.0 


13.0 


22.0 


2.0 


2.0 


7.0 


2. C 


3.0 


fc. " 


10.0 


31.3 


C.C 


0. 0 


6.0 


3.0 


3.0 


19.0 


17. r 


3 9.0 


4.C 


3.0 


13.0 


3.0 


1.0 


169. r 


12.0 


37.0 


1.0 


0.0 


12.0 


3.0 


3. 0 


102. 3 


11.0 


35. C 


2.0 


0 . 0 


11. r 


2.0 


2.0 


117.0 


9.0 


34.0 


l.C 


0.0 


10.0 


2.0 


3.0 


133.0 


13.0 


39.0 


13.0 


4. 0 


14. 0 


3.0 


7.0 


148.0 



NUMBER OB DATA PCINTS = 14 
RESULTS: 



ESTIMATE CF P 



81 

11.92C93 

27.11922 

2 . ^5082 
0. 89575 
6.64235 
2.22275 

3. '7*2 11 



P 2 

0.00456 
C. 06821 

0 • C 0 5 4 4 
0. CO 70 6 
0.03660 
0.00365 
-O.OC1 19 



ESTIMATE CF CCVARIANCE MATRIX SIGMA 



15.669 


9. 095 


4,728 


2.776 


3. 568 


3,183 


0 . «0 c 


9.095 


26.769 


-1.434 


0.710 


8.839 


6.3°6 


1.686 


4.728 


-1.434 


14.312 


3. 104 


0.379 


-0.975 


2. 02 p 


2.776 


0.710 


3.104 


2. 383 


0.515 


0.747 


0. 671 


3. 568 


8.839 


0.37 9 


0.515 


4.30 9 


2.357 


1. 551 


3.183 


6.396 


-0.975 


0,747 


2.357 


2.390 


1.151 


L. 80 5 


1.686 


2. 828 


0.671 


1.551 


1.151 


3. 237 



T test C p 


THE hYPG THESIS THAT A 


COMPONENT 


CF THE B2 VECTOR 


ECUALS 


ZERC AT A 


.l r SIGNIFICANCE 


LEVEL 


TEST STATISTIC 


RESULT 




C. 27445 


ACCEPT 




3. 14026 


REJECT 




C. 34273 


ACCEPT 




1.09012 


ACCEPT 




4.19994 


REJECT 




0.56178 


ACCEPT 




-0 .15716 


ACCEPT 
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ENGINE NUMBER: 516309 



CAT* : 



Y 1 


Y2 


Y 3 


Y 4 


7.-' 


31.0 


C.O 


C.O 


7 .' 1 


2 ).0 


8.0 


1.0 


6 . *•' 


1 9.0 


2.0 


0.0 


6 . 


21.0 


C.O 


0.0 


5 .r 


31.0 


C.o 


0 • 0 


5 • C 


26.0 


6 .o 


0.0 


1 5 . 0 


53.0 


5.0 


0.0 


3.'-' 


20.0 


8.0 


2.0 


A.r 


1 2.0 


C.O 


1.0 


NUMBER OF 


CAT A 


PCINTS = 


9 



Y 5 


Y 6 


Y7 


XI 


1C.C 


0.0 


3.0 


63. 0 


3.0 


o.c 


2.0 


32.0 


9.0 


0.0 


1.0 


26.0 


9.0 


• C.O 


O.C 


37.0 


7.0 


1.0 


1.0 


64. c 


9.C 


1.0 


2.0 


9C. 0 


17.0 


4 • C 


O.C 


ll c . o 


7.0 


0 . o 


?.c 


17.0 


4.0 


O.c 


0.0 


3 2.0 



RESULTS: 



ESTIMATE CF R 






82 



2 • 44 113 
P • 170 8 8 
3.< 32 8 C 
1. I I960 
4. 706 CC 
•I. 14859 
r. 74 821 



r- 



n 



• C 7 5 C 6 
. 34263 

noire 

■ 9 * C 1 2 6 6 
3 . n 78 4 3 
0.33404 
0.00472 



ESTIMATE CF CCVARIANCE m A jRIX SIG^A 



6.415 


13. 140 


1 • C Q 8 


-0.017 


5.261 


1.152 


- 0 . 


747 


13. 140 


49.284 


3.719 


0. 820 


13.528 


3.445 


-2. 


533 


1.098 


3.717 


14.2C6 


1 . 646 


2 . 3 °4 


0.795 


0. 


121 


-0. C17 


0. 8 20 


1.64 6 


"••,39 5 


-0.073 


0 . 1 8 C 


-c. 


2. '8 


5. 261 


13.529 


2,384 


- 9.073 


6.123 


n. qr p 


- 0 . 


4 82 


1. 152 


3.445 


0.7Q8 


A , 1 8 0 


9.90 8 


0.49 3 


_ r\ 

♦ 


63» 


-f . 747 


-2.533 


C . 1 2 1 


- 2?8 


-0.482 


-0.638 


1. 


4 0 C 



T TEST CF THE HYPOTHESIS THAT A 
COMPONENT CF THE 82 VECTCP ECUALS 



ZERO AT A .1C 



TEST STATISTIC 

2. 8 2 63 7 
4 . 6 c 78 1 
r . r 8989 
-1. c 2397 
3.C 2476 
4.62768 
0 . ? 3 T 3 6 



SIGNIFICANCE LEVFL 



RESULT 

REJECT 

REJECT 

ACCEPT 

REJECT 

REJECT 

REJECT 

ACCEPT 
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ENGINE NUMBER: 5 16094 



DATA: 



Y 1 


Y 2 


Y 3 


Y4 


Y5 


Y 6 


Y7 


XI 


c 

-- • ». 


16.0 


3. C 


9.0 


10. C 


l.C 


l.o 


le.o 


1C.C 


34.0 


7.0 


1.0 


11. f' 


2.0 


3.0 


94.0 


7.: 


2 B • 9 


6.0 


2. 3 


3 . 0 


l.C 


o.n 


62. C 


6,9 


2 3.0 


R.O 


2. 1 


B.~ 


2.0 


4.0 


93. 0 


l.i 


21.0 


l.C 


0.0 


8.0 


1.0 


l.C 


58. 0 


8. . 


32.0 


9.0 


0.0 


10.0 


1.0 


3.0 


88.9 


11.? 


3 1.0 


9.0 


0.0 


7.“' 


l.C 


2.0 


84. 0 


4.‘ 


27.0 


7.0 


1.0 


16.0 


2.0 


6 . C 


84. ^ 


7. 


24.0 


0.0 


0. 0 


5.0 


1.0 


1.0 


54.0 


MJNPLR OF 


LATA 


POINTS = 


c 











PESLLTS : 

ESTIMATE CF 8 
P 1 R 2 



4. 62-" 17 
12, 3°3 C 2 
-1.42774 
-r. 24627 
7.16616 
9. 49376 
- 1 • 0 ? c 1 1 



1 . C 3 6 6 3 
"'•27224 
7. C 99 12 
0 .U297 
0.C2866 
u.ll!9 A 
n. 64780 



ESTIVATE CF CCVARIANCc M ATR IX SIG^fi 



4.682 


3*213 


0.237 


-9. 817 


-4.107 


-0,551 


-2. 


«94 


3* 21C 


6.563 


1.813 


-0.012 


-1.492 


-0.429 


-2. 


173 


C.237 


1.813 


6.177 


0.322 


2.375 


-0.078 


C. 


539 


-C. 817 


-0.012 


''.32 2 


0. 737 


9. 257 


9.176 


- 0 . 


913 


-4. 1C7 


-1.492 


2.375 


r . 257 


19.486 


0. 803 


3. 


922 


-0.551 


-0.429 


-C. 078 


0.176 


0.89-3 


0.185 


0. 


459 


— 2 » 494 


-2.173 


0. 530 


- n . 0 1 3 


3. 922 


0.452 


2. 


373 



T TEST CF THE HYPOTHESIS THAT A 
COMPONENT C c THE B2 VECTCR ECUALS 
ZFRC AT A .10 SIGNIFICANCE LEVEL 



STATISTIC 


RESULT 


1.20353 


ACCEPT 


5. 57493 


REJECT 


2.81619 


REJECT 


1.L665C 


ACCEPT 


0.62926 


ACCEPT 


1 . c 5434 


REJECT 


2.19124 


REJECT 
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ENGINE NIJVRcr: 516949 



DATA : 



Y 1 


Y2 


Y 3 


Y4 


Y5 


Y 6 


Y7 


XI 


5* v 


22.0 


1.0 


1.0 


10. r 


2.C 


5 • C 


5 5.0 


7* ( 


19.0 


fi. 0 


2.0 


7.0 


2.C 


2.C 


8* 0 


12. ? 


2 1 . 0 


6 . r. 


C . o 


8.'' 


3 . 0 


2.0 


2 9. 0 


P." 


2 7.0 


7.0 


C . 0 


10 . c 


3.0 


2.0 


50.0 


6.F 


1 7.0 


2.0 


9. 0 


7.0 


i.r 


2.0 


5 5.0 


5* c 


24.0 


8.0 


2.0 


10.0 


2.0 


5.0 


8 5.9 


5.0 


19.0 


4 . C 


0.0 


6.0 


2.C 


3.0 


3 9.0 


6.C 


29.0 


0.0 


3.0 


13.0 


3 • C 


6.0 


86.0 


C r 

s # ** 


26.0 


3.0 


1.0 


11.9 


2.0 


4.0 


95.0 


8. > 


'3.0 


2.0 


1.0 


14.0 


1.0 


1 .C 


104.0 



NUMBEF CF TATA PCINTS = 10 



PFSILTS 



ESTIMATE CF P 



PI 

P. 14782 
14.?35ol 
6. P9968 
9. 36C23 
5*36915 
2.^7969 
2.C 7484 



B2 

-0.02389 
0.17268 
-0.C46 29 
0.91056 
0*06982 
-0.00627 
0.01857 



ESTIMATE CF CCVARIANCE W ATRIX S I G M A 



4. 891 


7.131 


U 96 0 


-0.600 


1.667 


°,499 


-1. 94? 


7. 131 


30.781 


1.478 


0.139 


6.50? 


0.216 


-5. 2 C 2 


C. 96C 


1.476 


7. 538 


?. 156 


-0.312 


0.29 7 


-9. 841 


-C.6CC 


o. 139 


r . ifo 


1.129 


0.69 7 


0.197 


1.036 


1.667 


6. 593 


-0.312 


0.697 


2.401 


0.526 


-0.062 


n,^.cc 


0.216 


0.297 


0.197 


0. e 26 


^ . 5 70 


C • 6 9 2 


-1.94 p 


-5. 292 


-0.841 


1.036 


-0.062 


0.60 2 


2. 82' 



T TEST OF THE HYPOTHESIS THAT A 



CrMDQNENT CF 
ZE R C AT A .10 



TFST STATISTIC 

-1.C085O 
2. c 0546 
-1. 57087 
C. 92769 
4.12975 
-C. 77438 
1.03131 



THE 02 VECTOR ECUAl S 
SIGNIFICANCE LEVFL 



RESULT 

ACCEPT 

REJECT 

ACCEPT 

ACCEPT 

REJECT 

ACCEPT 

ACCEPT 
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ENGINE NUMBER: 51675C 



DATA: 



Y 1 


Y 2 


Y 3 


Y4 


Y5 


Y 6 


Y7 


XI 


5.C 


11.0 


3.C 


0.0 


12.r 


2 • C 


2.0 


21.0 


6.C 


9.0 


2.0 


0.0 


12. C 


1.0 


2.C 


21.0 


6. 0 


13.0 


6.0 


1.0 


13.0 


3.0 


3.0 


27.0 


e.u 


28.0 


7.0 


1.0 


12. C 


1.0 


5.0 


61.0 


3.1- 


I 3.0 


2.0 


0.0 


4.C 


r* r> 

* * • W 


l.C 


96.9 


*.C 


16.0 


6.C 


1.0 


10.9 


1.0 


4.0 


13.0 


5.C 


15.0 


2.0 


1.0 


10.0 


c.c 


C.C 


15.0 


1C.C 


34.0 


9.0 


2.0 


12.0 


l.C 


3.C 


5*.Q 


8.0 


17.0 


6.0 


c.o 


12. n 


2 • G 


3.9 


22.0 



NUMBER OF DATA PCINTS = 9 



PESULTS: 

ESTIMATE CF B 
PI P2 



5.92P32 
12.5*119 
4. 3579* 
r. cc29^ 
13.C 3C24 
1 . 7*137 
2. *2961 



0.00514 

0.13476 

0.09868 

0.00207 

-0.06335 

-0.0146: 

0.00354 



F STI MATE CF CCVARIANCE MATRIX SIGMA 



5.536 


15.443 


4.30 1 


0.897 


4.700 


0.738 


1. 764 


15. * 43 


64.650 


16.711 


4.937 


11.656 


0. 20 8 


6. 840 


4.3C1 


16.711 


6.516 


1.272 


4.C23 


1.188 


3. 358 


0, 8 C 7 


4.937 


1.272 


0.5 68 


0.716 


-0.025 


0.376 


4. 70C 


11.656 


4.023 


9.716 


5.559 


1.384 


2. 324 


0,736 


0. 20 3 


1.188 


-0.025 


1.384 


C.°23 


0.736 


1.764 


6. 840 


3.358 


0. 376 


2.324 


0.736 


2. 594 



T TEST CF THE HYPOTHESIS THAT A 
COMPONENT CF THE B2 VECTOR ECUALS 
ZERO AT A .1C SIGNIFICANCE LEVEL 



TEST STATISTIC 


RESULT 


C. 1567C 


ACCEPT 


1.20225 


ACCEPT 


9.24397 


ACCEPT 


C. 19728 


ACCEPT 


-1.92719 


REJECT 


-l.C 9022 


ACCEPT 


f. 15775 


ACCEPT 
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ENGINE NUMPER: 51657* 



DATA: 



Y 1 


Y2 


Y 3 


Y4 


Y5 


Y 6 


Y7 


XI 


7 • i 


36. C 


r r* 
v • ' 


0. 3 


9.0 


1.0 


6.0 


2 2.o 


5.( 


37. 0 


n 

w# • < 


0 . ') 


I2. r 


2.C 


2.0 


64.0 


6. ~ 


24. ) 


3.0 


c.o 


11.0 


2 . C 


2.0 


5 6.0 


17. ) 


29.0 


4.0 


c, 0 


9." 


• 2.0 


1.0 


28.0 


16.’ 


3 5.0 


9.0 


1.0 


9.0 


c r> 


4.C 


36. 7 


1 6 . 


51.0 


6. 9 


2. 0 


17." 


4.0 


4.0 


112.7 


10.0 


3 9.0 


8.0 


1. 1 


i r\ -o 
i 


4 . C 


4.C 


8 5.0 


3 3.0 


3 4.0 


2.0 


0. > 


9.0 


5.0 


3.0 


63.o 


C t 

~ • 


25. 0 


7. C 


1.0 


L r 

• V 


0.0 


-1.0 


27.0 



MJMRFP OF TATA PCINTS = 9 



RESULTS : 

ESTIMATE cf b 

Ml B2 

8 . 7 8 7 6 C 1.03202 

24.7697F '>.19532 

3.14465 0.01951 

-(.. 16 926 G. 013 13 

5. •.*1866 j. 09021 

1.35( 9 8 0.03127 

2.18423 1.C1276 



ESTIMATE CF CCVARIANCE MATRIX SIGMA 



27. £19 


3.221 


6.475 


1.022 


1.557 


5.828 


1. 


510 


3. 220 


16. 5 14 


-3.422 


0.356 


3.89e 


0.78 5 


6. 


<+78 


6.475 


-3.422 


11.253 


1.717 


-2.551 


1.427 


-1. 


0 85 


1. 022 


0.356 


1.717 


0.423 


-0.093 


0.016 


o. 


04 7 


1 . 557 


3.898 


-2.551 


-o.093 


2.936 


-0. SCI 


C. 


22? 


5. 82 8 


0.785 


1.427 


0 .016 


-0.801 


2.630 


1. 


123 


1.510 


6.478 


-1.085 


0. r 47 


0.227 


1.123 


3. 


671 



T TEST CF THE HYPOTHESIS THAT A 
COMPONENT CF THE R2 VECTOR ECUALS 



ZFRC AT A .1C 



TFST STATISTIC 

f . 5 1 8 8 1 
4.13930 
C .49718 
1.72416 
4.49939 
1.64798 
C. 5S922 



SIGNIFICANCE LEVEL 



RESULT 

ACCEPT 

REJECT 

ACCFPT 

ACCEPT 

REJECT 

ACCEPT 

ACCEPT 
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ENGINE NIJMRER: 516534 



DATA : 



Y 1 


Y2 


Y3 


Y4 


Y5 


Y 6 


Y7 


XI 


11.0 


27.0 


5.0 


0.0 


8.0 


1.0 


4.9 


17.0 


14. u 


32.3 


7.0 


1.0 


14.0 


1.0 


3.0 


19. " 


15.0 


4 3.0 


6.0 


0. 9 


17.0 


2.0 


4. r 


8C.0 


15.0 


42.0 


7.0 


2.0 


15.0 


3 • c 


1.6 


84.0 


14. n 


28.0 


7.0 


3.0 


8.0 


2.0 


3.C 


24.0 


18. c 


27.0 


1.0 


3.9 


7.0 


2.0 


2.C 


24.0 


14.0 


34.0 


10.0 


1.0 


9.0 


2.0 


6.0 


9.0 


17.0 


3 9.0 


1.0 


0.3 


15.0 


2.0 


5.0 


54.0 


14. T 


38.9 


8.0 


0. 3 


10. C 


2.0 


3.0 


32.0 


NUMBER OF 


DATA 


FCINTS = 


9 











FFSULTS: 

ESTIMATE CF e 



B1 

13.79272 
27. 26 °3 5 
6.4 l*+ c 9 
1.36354 
7.34557 
1.366 15 
4. 1 34P7 



R2 

0.02293 
0.18827 
-0.0167? 
-0.00662 
0. 10755 
0.01387 
-O.C1812 



ESTIMATE CF CCVARIANCE MATRIX SIGMA 



< « 1 C 5 -0.63 8 

0.638 14.603 

3.898 6.493 

1.182 -3.101 

r.42f 5. €78 

C ♦ 385 9.176 

0.298 3.199 



-3.898 1.132 

6.490 -3.101 

10.260 -9,495 

-C.495 1.802 

1.293 -1.718 

0.317 0.383 

0.430 -1.313 



-0.425 3.385 

5.078 9.176 

1.293 0.317 

-1.718 0.383 

6*061 -0.688 

- 0.688 0.242 

1.617 -0.142 



-C* 298 
3. 199 
C.43C 
-1.313 
1.617 
-0. 142 
2. 312 



T TEST OF THF HYPOTHESIS THAT A 
COMPONENT CF THE 82 VECTCR ECUALS 
ZERO AT A .10 SIGNIFICANCE LEVEL 



TF ST STATISTIC 


RESULT 


C. 89167 


ACCEPT 


3.88142 


OBJECT 


-€. 41124 


ACCEPT 


-0.38869 


ACC C PT 


3.44178 


REJECT 


2. 22170 


REJECT 


-C. 93863 


ACCEPT 
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ENGINE NUMBER 



5) 90 9.1 



DATA: 



Y 1 


Y 2 


Y 3 


Y4 


Y5 


Y 6 


Y7 


XI 


4 , - ■ 


17.0 


6.0 


1. 7 


5.o 


1.0 


1.0 


57.0 


1 °.r 


2 9.0 


1.0 


0.0 


1.0 


3 • C 


1.0 


91.0 


2. ' 


19.0 


0.0 


1.0 


4.o 


0.0 


0.0 


63.0 


7.o 


3 2.0 


5.0 


0.0 


9.0 


2.0 


o. r. 


180. 0 


P.L 


39,0 


0.0 


o. 1 


9.0 


3.0 


2.0 


209. 0 


r r 
- • ' 


21.0 


0.0 


0.0 


5.0 


o.o 


l.o 


237.0 


7.0 


27.0 


4.0 


c.o 


5.0 


3.0 


1.0 


148.0 


3.<" 


26. 0 


6.0 


0.0 


l. r 


l.C 


l.C 


110.0 


6 . <. 


21.0 


8.0 


0.0 


6.0 


1.0 


0.0 


119.0 


9.o 


24.0 


2.0 


0.0 


6." 


2.° 


2.0 


106. 0 



NUMBER OF DATA FCJNTS = 10 



RESULT?: 



ESTIMATE CF fi 
PI B2 



7.O20C2 
17. 745< A 
3.2^251 
f.775Bf 
3.i 582C 
1.2371E 
C.f 56( C 



0.35879 
0.0065C 
-C .00^37 
0.01973 
0 . 0 02 7 5 
U. 00261 



ESTIMATE CF CCVARIANCE MATRIX SIGMA 



25. 


50 C 


15.287 


-3.121 


-1.003 


-5.738 


4.377 


1. 


127 


1 5 • 


287 


35.331 


7.985 


-0.818 


3.663 


5.959 


1. 


431 


— *3 

— • 


121 


7 • c 8 9 


11.689 


-0. 158 


5.135 


1.226 


0 • 


193 


-1. 


CO 3 


-0.813 


-0.158 


0. 122 


0.0 5 5 


-0,226 


-0. 


053 


_ e 

* 


738 


3. 663 


5.135 


0.C5 5 


4.661 


0 . 00 2 


o. 


O^l 


^ * 


377 


5.959 


1.226 


-0 .226 


0.00 2 


1.519 


o. 


420 


l. 


127 


1.431 


0.193 


-C.053 


O.OC! 


0. 4 20 


r'. 

‘ • 


584 



T TEST OF THE HYPOTHESIS THAT A 
COMPONENT CF THE P2 VECTOR FCUALS 



Z F R C AT A .10 



TEST STATISTIC 

-r .C 1546 
1.79515 
C. 34506 
-2. 27294 
1.65826 
0.41507 
C. 61910 



SIGNIFICANCE LEVEL 



R C SULT 

ACCEPT 

ACCEPT 

ACCEPT 

REJECT 

ACCEPT 

ACCEPT 

ACCEPT 
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ENGINF NUMBER: 516592 



CAT A : 



Y1 


Y2 


Y 3 


Y4 


Y 5 


Y 6 


Y7 


XI 


5.r 


21.0 


U • 


0.0 


12.0 


2.0 


4.C 


89.0 


3.? 


11.0 


5 . C 


1 . 0 


6.0 


2 • C 


2. C 


30.0 


5 • i 


2 2.0 


4." 


0.0 


10.0 


4.0 


3.0 


59.0 


5.< 


19.0 


5.0 


1.0 


9.* 


3.0 


2.0 


83. ° 


6. ° 


17.0 


2.C 


1.0 


10. o 


2.0 


4.0 


95.0 


9.0 


18.0 


c.c 


0 . 0 


11. n 


3.0 


5.0 


63.0 


7. 0 


15.0 


6.0 


1.0 


8 .* 


2.0 


4.0 


37.0 


6. i 


16.0 


2.0 


1.0 


10. c 


2.0 


c.o 


33.0 


6.0 


18.0 


4.0 


2. 0 


4.0 


1.0 


o.C 


34.0 



NUMBER OF TATA PCINTS = 9 

RESULTS: 

ESTIMATE CF B 



01 


R2 


4.66621 


0.C2678 


13.2T392 


0.07257 


5.39C8C 


- 0 • 0 3 c 2 3 


1.4-3C99 


-C.O 11 24 


5.22276 


0.063C9 


1.73505 


o.O If 30 


r. '■••383 2 


0- . C 3 8 3 5 



ESTIMATE CF CCVARI ANC6 MATRIX SIGMA 



3.66 8 


1.792 


-2.649 


-0.418 


1.87C 


- 0 . 3 r 8 


1. 302 


1. 792 


S.204 


- 0! . 9 9 5 


-0.525 


1.222 


•''.944 


-0. 828 


-2. 64 9 


-0.995 


4.367 


0. 548 


-2.358 


r . 264 


-0. 362 


-r, 4ie 


-0.525 


0.648 


0,410 


-1.055 


-0. 387 


-C. 619 


1. F7C 


1.222 


-2.358 


-1,055 


4.196 


0.83 2 


1.227 


-G. 306 


0. 944 


0.264 


-0.387 


0.832 


0.775 


0.409 


1.3C2 


-0. 623 


-0.362 


-0.619 


1.227 


0.409 


2. 579 



T TEST CF 


THE HYPOTHESIS THAT A 


c cmfonent 


CF THE P 2 VECTOR 


EQUALS 


ZFRC AT A 


.10 SIGNIFICANCE 


LEVEL 


TEST STATISTIC 


RESULT 




1.5 2799 


ACCEPT 




1. 87332 


ACCEPT 




-1.33041 


ACCEPT 




-1. 29033 


ACCEPT 




2.26479 


REJECT 




r . 85980 


ACCEPT 




1.75595 


ACCEPT 
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ENGINE NUMBER : 515931 



DATA : 



Y 1 


Y2 


Y3 


Y4 


18. O 


35. G 


0.9 


1.0 


10.0 


27.0 


e.c 


o.o 


6.0 


34. ) 


3.0 


0.0 


1 c . r, 


37. D 


7.0 


1. 3 


P.v 


2 Q 


5.0 


0. 0 


12.f 


34.6 


6. 0 


1.0 


5.0 


23.0 


5.0 


0.0 


lo. r 


33.0 


2.o 


l.o 


6.6 


27.0 


o . ) 


0. 0 


5. r 


21.0 


4 • 0 


0. o 



NUMBER OF DATA FCIMS = 10 



Y 5 


Y 6 


Y7 


XI 


I l.o 


2.0 


3.0 


60. 0 


10.0 


3.C 


2.0 


36. 0 


11.0 


2.0 


3.0 


31.0 


14.0 


• 2.0 


3.C 


120. C 


10.0 


2.0 


2.9 


09.0 


ll.o 


o.Q 


5.0 


89.9 


8.0 


e.c 


l.o 


18.0 


1 1 . 0 


0.0 


1.0 


75. C 


8.0 


0.0 


2.0 


51. 0 


7.0 


c . 0 


l.C 


9.0 



F FSLLTS: 



B 1 



ESTIMATE cf p 
B2 



6.43751 
23.349^7 
l. c c-7ei 
-C.184?E 
7. 4 8260 
C • 7 5 6 C 6 
1.29668 



3.04747 
,.11599 
0 .029^8 
0.01004 
0. n 4496 
C . 0 0 5 c 1 
0.01720 



ESTIMATE CF COVARIANCE MATRIX SIG W A 



13. 


719 


8.70 4 


-4 . Q04 


). Q6 6 


2.165 


1.658 


1. 


07 c 


e. 


704 


13.568 


-2.737 


G. 655 


3.884 


1.755 


2. 


22 5 


-4. 


CP4 


-2.737 


6.4PP 


-o. 144 


-0. 1 5 r 


-0.896 


0 . 


20 8 


C, 


566 


0.655 


-o. 144 


<'.155 


0.178 


-0. 135 


9. 


102 


2 * 


165 


3.8 84 


-0.150 


9. 178 


1.714 


n .9fi2 


' • 


47 6 


1. 


698 


1.755 


396 


-9.135 


3.9R2 


1.56 2 


0 . 


1°2 


1. 


879 


2.225 


C . 2 0 8 


3. 10 2 


0.47 3 


0. 192 


1. 


338 



T TEST CF THE HYPOTHESIS THAT A 
COMPONENT CF THE P2 VECTCP ECUALS 
7FRC AT A .10 SIGNIFICANCE LEVEL 



TEST STATISTIC RESULT 

1.37266 ACCEPT 
3.372Q3 REJECT 
1.22274 ACCEPT 
2.72805 REJECT 
3.67812 REJECT 
0.60639 ACCEPT 
1.59321 ACCEPT 
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ENGINE NUMBER: 5 1 5 6 0 e 



DATA : 



Y 1 


Y 2 


Y 3 


Y4 


6 . C 


33.0 


4 . C 


1.0 


8 .* 


4 >. 0 


1.0 


0.0 


4.0 


?l.o 


10 . c 


0.0 


7 • •" 


31.0 


7. r 


0 . 0 


K • 0 


37.0 


l.v 


1.9 


2 . f 


3 3.0 


3. C 


1.3 


4.f 


21.0 


3.0 


2.0 


7 • ' 


27. <> 


4.0 


0 . 1 


6 .. 


3 1.0 


9.0 


1.9 


NUMPLR OF 


DATA 


POINTS = 


9 



RESULTS : 

FSTIMATF C p p 



81 


R2 


fc . 2 c 5 8 2 


-”’•00061 


3( . 1 34 7 C 


0 . *’-2774 


8 • * 4 7 


-C.C611C 


r. 3 5927 


0.w0557 


11.71527 


0.C7357 


G. 78598 


0. Gil 92 


1.55565 


0.C2817 



Y5 


Y 6 


Y7 


XI 


16.0 


1,3 


3.0 


2 2.3 


22 . n 


2 


4.0 


66.0 


18.3 


l.C 


4.0 


58. r 


14. 0 


1.0 


2.0 


2 8 . 0 


21.0 


2.0 


6 . C 


112 .C 


15.0 


2 .o 


3.0 


8 5.0 


12 . C 


l.C 


0.0 


63. 0 


12 . C 


3.C 


5.0 


36. 


12 . C 


0 • ♦.* 


l.C 


27.0 



ESTIMATE CF CCVARIANCE MATRIX SIGV4 



7.079 
1C, 69* 
-3.229 
-0.473 
4 . 967 
C.452 
2.7CC 



10.684 

39.1Q7 
-4.374 
-2.018 
18.775 
0.274 
6 . 375 



-3.229 
-4.374 
8 . 3*30 
-C.359 
- 1 . 11 ° 
-1.335 
-0.585 



—*'.473 
-2.C 18 
-0,350 
0. 539 
-1.380 
-3. 30 7 
-0.973 



4.967 
13.775 
-1.119 
- 1. 330 
11.145 
0. r 6 c 
3.155 



0.452 2.7GC 

0.274 6.375 

-1.335 - C . 585 

-0,307 -C.973 

C.069 3.155 

0. 741 1.015 

1, ~15 3.295 



T TEST OF THF HYPCTHFSIS THAT A 
OPPONENT CF TEE 82 VECTOR EQUALS 



ZFRC AT A .1C 



TEST STATISTIC 

-C.C1954 
0.37673 
-1.6C28C 
f. 64801 
1. 83347 
1. 18436 
1.32551 



SIGNI c ICANCF LEVEL 



RESULT 

ACCEPT 

ACCEPT 

ACCEPT 

ACCEPT 

ACCEPT 

ACCEPT 

ACCEPT 
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ENGINE NUMBER: 51638C 



rATA : 



Y 1 


Y 2 


Y3 


Y4 


Y 5 


Y6 


Y7 


XI 


5.0 


19.0 


6.0 


C. 9 


9.0 


1.0 


4.0 


60.0 


U . 0 


3 5.0 


3.0 


1.0 


15. C 


l.C 


1.0 


58.0 


18.0 


4 3 . 0 


5.0 


2.0 


13.0 


2.0 


5.0 


5 2.0 


4.»* 


24.0 


7.0 


2.0 


10. 0 


l.C 


l.C 


12.0 


6 • 0 


27.0 


9.C 


l.o 


10.0 


l.C 


l.C 


29.0 


6. ( 


25.0 


1.0 


1. J 


12.0 


l.C 


l.C 


41.0 


5.0 


26.0 


6.0 


0.0 


12.0 


1.0 


2.0 


6fc.o 


E r. 

- • 


19.0 


7 . C 


3.o 


8.0 


l.C. 


5.0 


6.0 


*. ) 


16. G 


5.0 


o.O 


4 • ° 


1.0 


0.0 


39.0 


5.0 


1 5 . 0 


5 • C 


1.0 


5.0 


l.C 


n .c 


31.0 


3.-: 


19.0 


2. A 


0.0 


6 • 


1.0 


2.0 


18.0 



NUMBER OF TATA PCI NTS = 11 



RESULTS: 

FSTIMATE C R B 



PI 

3.26027 
17.740C 8 
5- 9444? 
1.97730 
6 . 2*200 
C.958C6 
1.P7213 



P2 

0.08528 
0.17684 
-0.02279 
- A .C2609 
0.08577 
0.00355 
0.CP34 1 



FSTIMATE CF COVARIANCE MATRIX SIGMA 



16. 


545 


29.482 


C.0 57 


2. 681 


6.526 


1. 145 


3. 756 


28. 


482 


65. 145 


0.463 


4. 880 


20.665 


1.7e5 


5. 392 


C. 


r 57 


0.463 


5.864 


0.507 


-0.049 


0.027 


0.924 


2. 


681 


4. 880 


0.50 7 


0.801 


l. Q 09 


0.153 


0. 929 


6. 


526 


20.665 


-C. 049 


1.909 


10.062 


0.255 


1.644 


1. 


145 


1.785 


G.02 7 


0. 153 


0.255 


0.095 


C. 328 


3. 


75 6 


5. 392 


C. R 24 


a. c 2 9 


1.644 


C.328 


3.772 



T TEST OF THE HYPOTHESIS THAT A 
COMPONENT OF THE R2 VECTOR ECUALS 



ZERC AT A .10 



TEST STATISTIC 

1.34267 
1.40306 
- 0.60264 
-1. 86709 
1.73153 
0.73587 
0 .1125 6 



SIGNIFICANCE LEVEL 



RESULT 

ACCEPT 

ACCEPT 

ACCEPT 

REJECT 

ACCEPT 

ACCEPT 

ACCEPT 
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ENGINE NUMBER: 515515 

fATA : 



Y! 


Y 2 


Y 3 


Y4 


Y 5 


Y 6 


Y 7 


XI 


IF... 


4 1 . 0 


7 . 0 


C. 1 


13.0 


2.0 


6.0 


42. 9 


1 8 . o 


4 5.9 


1.0- 


9 . 9 


11.0 


3 • C 


5.0 


55. C 


17. ' 


A 5. ) 


3.9 


0. 9 


11.9 


2.0 


4.G 


5 8. 0 


12.0 


9 6.1) 


9. ^ 


1.0 


12." 


2 . C 


4 . C 


37.9 


1 2 . ( 


37. ) 


7. o< 


5.9 


1 0 . 0 


2.: 


l.C 


24.0 


1C.. 


3 ?, <•* 


3.*' 


1.9 


8.0 


2.0 


0.0 


24.0 


c . ' 


34,6 


1.9 


1. J 


13.9 


1.9 


3.0 


62. 9 


18.! 


a 3 . ') 


2.1 


0.9 


1 5 • C 


2.0 


5 . C 


76.9 


11. 


3 6.9 


6.0 


9. 3 


8 . 0 


2.0 


2.9 


3.0 



NUMBER HF TATA priMS = 9 

RESULTS : 

ESTIMATE CF P 



Q 1 


82 


If. 12277 


0.08696 


32. 5r 766 


3,16124 


7 • 6 C 92 5 


- 0 . 07738 


r . - 577 c f 


— r # 9~o 9 a 


7.577F1 


C.C8699 


2. . 79oo; 


-0.9.-166 


C. c 7r 21 


9 . 0 5 5 R 2 



FSTIMATF CF CCVAPIANCF MATRIX sigma 



11.693 
14. r it 

2. 153 

-1.A73 

C. 46 c 
1.375 

3. 9^4 



14, r 66 
20.693 
1.44 J 
-2.194 
■ 9 . 05 7 
1.733 

4. "62 



2. 153 

1.449 
6 . a 2 9 
-0.044 
2.591 
-0.077 
2. 1 A 7 



- 1 . A 7 3 
-2.194 
- 9.944 

9. 285 
-0.946 
. 144 
-0. 397 



'■.465 

0. 957 
2 • 5 Q 1 

- ? • r 4 6 
1.893 
-0.2CC 

1. 182 



1.375 
1. 733 
- 0.977 
-9.144 

-0.2GC 
0.2 84 
9,342 



3. 944 
4.96 2 

2. 14 
-C.39 
1. 18 
0.34 
2.71 



T TEST CF THE HYPOTHESIS THAT A 
COMPONENT CF THE 82 VECTCR ECIJALS 
7FRC AT A .10 SIGNIFICANCE LEVEL 



TEST STATISTIC 


RESULT 


1.6 3 180 


ACCEPT 


2. 30081 


PE IFCT 


- 1. c 717? 


REJECT 


-0. 11594 


ACCEPT 


a. 14212 


REJECT 


-C. 2 J33C* 


ACCEPT 


2. 13941 


PE JECT 
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ENGINE N U M R E R : 516731 



CATA : 



Y 1 


Y2 


Y 3 


Y 4 


Y5 


Y 6 


Y7 


XI 


e.c 


16.0 


e.r 


I."' 


14.o 


2.G 


5 • C 


54. (, 


ll.f 


29.0 


6.c> 


2.0 


16.0 


3.0 


6 • (; 


8 3.0 


13. C 


21.0 


B.C 


2. 0 


9.0 


l.C 


3.0 


57.C 


6.v 


2 0.0 


6.0 


3. ? 


9.0 


O.C 


5 • C 


3 3.0 


C 

^ • ‘ 


20.0 


8.0 


0.0 


8.0 


1.0 


1.0 


8 7. A 


10.4 


22. 0 


5.0 


1.0 


13.0 


2.0 


3.0 


27.0 


13. C 


24.0 


5.0 


1.0 


7.0 


1.0 


0.0 


2 2.0 


5. 


2 3.0 


7.0 


2.0 


9.0 


2 • C 


4.0 


99. 0 


5.0 


15.0 


5.0 


1.0 


8.0 


l.C 


O.C 


31.C 



NUMBER PF DATA PCINTS = 9 

RESULTS: 

FSTIMAT6 CF P 



81 


P 2 


10.39656 


-C. 03564 


ie.452P< 


7.04852 


4.6C35C 


O.C2 Q 55 


1.45355 


-o. con i7 


c . 1 P P 3 4 


0.C209C 


0.67471 


C . * l 4 0 5 


1.45051 


C. 02829 



FSTIMATE CF 


CCVAP I ANCE M AT R I X SIGMA 






11.962 


9. 873 


C.45C 


'). 312 


3.52 I 


1.224 


1. 248 


9. E73 


17.879 


-1.115 


1.230 


3.9 84 


1.429 


2. 261 


0 . 4-50 


-1.115 


0.817 


0.021 


-1.399 


-0.523 


-0. 369 


C.312 


1.237 


'->.021 


0 .889 


0.384 


- 0. 109 


1. 433 


3. 521 


3.984 


-1. 399 


0. 384 


11.011 


2.24-3 


5.436 


1. 224 


1.429 


-0.523 


-0 • 109 


2.243 


0.700 


0.621 


1. 24 8 


2.261 


-l.369 


1.433 


5.436 


n .621 


4.951 



T TEST GF THF HYPOTHESIS THAT A 
COMPONENT CF THE R2 VECTOR EQUALS 
Z r R C AT A .1C SIGNIFICANCE LEVEL 



STATI STIC 


RESULT 


P. 84223 


ACCEPT 


C. c 3905 


ACCEPT 


2. 67191 


REJECT 


C. r 1442 


ACCEPT 


0.51491 


ACCEPT 


1.37245 


ACCEPT 


l.C 3920 


ACCEPT 
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ENGINE NIJVEER: 515447 



DATA : 



VI 


V2 


V 3 


Y4 


Y5 


Y 6 


Y7 


XI 


10 . V 


32. 5 


9. • 


o.o 


11.'' 


l.c 


2.0 


65. ~ 


12.. 


M . ! 


2.'« 


3. n 


1 5 . r 


2. 


2M 


82.0 


P. 


26. ^ 


4 . • » 


0. 1 


9.0 


6.0 


O.C 


46. 


7. '' 


2 8.-' 


4 • ' 


1.0 


11.9 


l.r 


2. C 


5 4 . 0 


P.o 


2 6.0 


l.c 


7. o 


10. 


1.0 


?. r 


6 : . c 


5.'- 


27.0 


3.o 


1. 1 


9.0 


l . r 


0 . 0 


7 1.0 


4 . 


21.9 


4.0 


*>.0 


8.2 


1.0 


C . 0 


2 6. o 


6. : 


27. •) 


7. ~ 


1.0, 


9.-' 


l.c 


3.0 


2 6. 0 


12. '• 


39.0 


3. r 


r 

v. • . 


12.' 


fC 

• * » 


2 . c 


8 2.0, 


12.1 


26.0 


2.0 


o.o 


8. ^ 


4 . C 


'.C 


70.-' 



MJNPEP OP TATA PC I MS = l n 



RESULTS : 

1ST! MATE CP P 
PI P 2 



1.992C1 
It.? 5591 
6.64628 
-O.l 0295 
6. 02 c 82 
C. P4f c 4 

1 • Z 5 62 8 



i ioio 

0.22241 

-0.04719 

. . 013PC 
■'•06485 
' • vl 2 5 " 7 
0.00247 



ESTIMATE CF CCVARIANCF MATRIX SIGMA 



4.4(4 


3.621 


0.845 


-0. 357 


0.825 


2.06 5 


0.673 


3.621 


2 1.914 


4.162 


1.496 


7.358 


", 20 3 


3. 519 


C. P* 5 


4.162 


5. 576 


— 0. 847 


C.625 


-0.787 


0. 729 


-o. 357 


1.496 


-0.847 


1.111 


1.133 


-1.086 


9 . 799 


C. P25 


7.35 8 


0 • P 2 5< 


1.133 


3,242 


-0.832 


1.460 


2. 0 65 


G. 2C3 


-<'.7P7 


-1.086 


-0.832 


3,970 


-1.054 


0.67? 


3.519 


-.729 


o. 7 ’1Q 


1.4 5C 


-1.954 


1. 797 



T TEST CF THE PYCCTHESIS THAT A 
CC W PCNENT CF THE P2 VECTOR ECU ALS 



ZERC AT A .1C 



TEST STATISTIC 

3 • 2 J ? 2 7 
2. 699-91 
-1. 21977 
C .79850 
2.19847 
C . 76801 
0.1124? 



SIGNIFICANCE level 



RESULT 

REJFCT 

REJECT 

ACCEPT 

ACCEPT 

REJECT 

ACCEPT 

ACCEPT 
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ENGINE NUMBER: 516048 



TATA : 



Y 1 


Y 2 


Y 3 


Y4 


6 . ■ 


24.0 


1.0 


0.0 


c . ( 


/ *6. 0 


4.^ 


2.0 


A • r 


24.0 


6 . v 


1.0 


9.0 


4 6. 0 


5.C 


1.0 


L.i 


19.0 


3.0 


2.0 


6. ‘ 


39.0 


6.0 


l.Ci 


A .f’ 


29.0 


3.0 


2.0 


1C.' 


c 2 • 0 


8.0 


1.0 


2 . • ■ 


(\J 

O 

• 

c 


C.C 


0. 0 



MJNPEP OF TATA POINTS = 5 



Y 5 


Y 6 


Y7 


XI 


P.o 


0.0 


2.0 


30. t 


12.0 


2 . o 


3 • C 


120.0 


6.o 


o.C 


0.0 


90. " 


13.''' 


2.0 


4.0 


96.0 


5.0 


1.0 


1.0 


55.0 


11.0 


l.C 


1.0 


55.0 


9.C 


1.0 


3.0 


3 3. ° 


16.C 


1.9 


2.0 


96.0 


8.0 


l.C 


0.0 


2 °.n 



FESULTS: 

FSTIMATE CF B 
R 1 P2 



1. FA All 
1 5.05L 2< 
C .53192 
•' . A833C 
6. ! FI 3 3 
f . 29463 
C. 83647 



0 • 0 6 0 r A 
0.22459 
C.C5134 
0.C0929 
C.C53C 9 
9 • C 1 r 4 4 
0 . C 1 3 9 3 



ESTIMATE CF CCVARIANCE MATRIX SIGMA 



4.123 


14. 290 


0.552 


-r . 304 


5.385 


0.319 


1. 901 


14. 29C 


92. 544 


5. <<28 


-1.773 


30.5^6 


3.206 


6. 096 


C • 9 5 2 


5.428 


3.96 7 


0.088 


2.277 


-C.418 


- 0 . 366 


-C . 3C-A 


-1.773 


0 . C 8 8 


0. 585 


-C.616 


r >. 159 


C. 290 


5.385 


30. 5 A- 6 


2.277 


-0.616 


10. 52C 


9, 843 


2. 251 


C. 315 


3.206 


-0.418 


0. 158 


0.843 


0.42 8 


C. 523 


1. 9C1 


6 . C 9 6 


-C.36S 


C. 290 


2.251 


0.523 


1. 967 


T 


TEST CF 


THE HYPOTHESIS THAT A 






CC 


*P0NENT 


CF THE 82 


VECTOR 


ECUALS 







ZFRC AT A .1° SIGNIFICANCE LEVFl 



TEST STATISTIC RESULT 

2.83524 REJECT 
2.23854 REJECT 
2 • 4 7 1 46 REJECT 
1.15508 ACC E °T 
1.56545 ACCEPT 
1.52994 ACCEPT 
C. c 5259 ACCEPT 
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ENGINE NUMBER: 516063 



CATA : 



Y 1 


Y2 


Y 3 


Y4 


Y5 


Y 6 


Y7 


XI 


2.;' 


11.0 


3.C 


C.O 


6.0 


6 . 0 


2.0 


C.o 


7. r 


2 0.0 


0 • 0 


0.0 


12.0 


1.0 


4.0 


71.0 


6. r 


18.0 


2.0 


0.0 


10.0 


1.0 


2.0 


34. C 


6. 7 


19.0 


C.O 


1.0 


°.o 


1.0 


4.o 


40. " 


2. <*• 


1 3.0 


8.0 


1. ) 


8.0 


0.0 


2.0 


27.0 


9. r ‘ 


23.0 


l.C 


0.0 


14. Of 


1.0 


2.o 


0.0 


4.0 


1 7.0 


9.0 


1.0 


9.0 


2.0 


4.0 


24.0 


5.C 


20.0 


5.0 


1. ) 


9.0 


7.0 


2.0 


17.0 


7.r 


32.0 


2.0 


C.O 


14.0 


7.0 


4.0 


65. C 



NUMBER CF TATA PCINTS = 9 

R f SILT $ : 

FSTIMATE CF P 



Q 1 


B2 


6.53252 


0 . 0391 ? 


15.2696* 


0.12472 


4. 63345 


-0.04209 


f . 53959 


- 0 • 0 0 3 r 9 


E. 66085 


C. 04655 


3. 1 1C6 C 


-0.07718 


1 • C 4 C e«. 


G.C 3 n 4C 



ESTIMATE CF COVARIANCE MATRIX SIG^A 



6. 


387 


9. 302 


-5.951 


-0.387 


5.089 


-1.321 


0. 


356 


9. 


3C2 


30.475 


-5. 324 


-0 . 850 


12.190 


6.961 


c • 


450 


-5 « 


55 ! 


-5.324 


11.299 


1. 145 


-3. 342 


-0.027 


c. 


253 


-C. 


387 


-0.850 


1.145 


0. 311 


-0.673 


-0.238 


o. 


131 


• 


089 


12. 19’^ 


-3.342 


-0.673 


6.829 


-0.17C 


-0. 


CIO 


-1. 


321 


6.961 


-0 .027 


-0.238 


-0.170 


9.519 




00 2 


C . 


356 


0.45 C 


0.253 


0.131 


-0.0 10 


-0.0C2 


o. 


60 6 



T TEST CF 


THE HYPOTHFSIS THAT A 


CCVPONENT 


CF THE e2 VECTGR 


ECUALS 


ZERC AT A 


.1C SIGNIFICANCE 


IEVFL 


TEST STATISTIC 


RESULT 




C. c 2937 


ACCEPT 




1.6)218 


ACCEPT 




-C. 83798 


ACCEPT 




-C. 39369 


ACCEPT 




1.27409 


ACCEPT 




-C. 16505 


ACCEPT 




2.76992 


REJECT 
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ENGINE NUMBER: 539064, 

DATA : 



Y 1 


Y 2 


Y 3 


Y4 


Y5 


Y 6 


Y7 


XI 


14.0 


45. 1 


l.n 


1.0 


22.0 


2.0 


4.0 


240.0 


11.0 


31.0 


2.0 


1. ) 


11.0 


3 • C 


5.0 


103.0 


9.1 


36.0 


6.0 


0.0 


13.0 


4.0 ' 


2.0 


90.0 


10.0 


3 5.0 


7.0 


1.0 


11.0 


4.0 


4.0 


67.0 


1C. o 


24.3 


7.C 


1. 0 


3.0 


5 . C 


1.0 


30. C 


1 4 . ’ : 


42.0 


8.0 


2.0 


19. r 


2.0 


6.0 


210.0 


14. 


4 3.0 


e.o 


3.0 


19. C 


2.0 


6.0 


217.0 


ll.o 


3 5.0 


9.0 


0.0 


13.0 


2.0 


5.0 


150.0 


e. i 


31.0 


5 • c 


1.0 


10 . c 


0 . 0 


3.0 


60.0 



NUNBfcP OF DATA POINTS = 9 



RESULTS: 

ESTIMATE CF R 
B 1 R 2 



7.72752 
25.42332 
6. 27336 

f . 364 84 
6.15189 
3.73C31 
1.88632 



0.02655 
0.07985 
-0.00,297 
C. 00576 
0.06053 
-C. 009 20 
0. Cl 63 3 



F ST I MAT t CF CCVARIANCE M AT R I X SIGMA 



c. 


QC( 


-1. 163 


C.299 


0.483 


-0.108 


0.886 


0. 


151 


-1. 


163 


8.229 


0.152 


-0. 098 


2.517 


-0.063 


- 0 . 


150 


0. 


2 90 


0.152 


8.635 


0.541 


-1.064 


0.358 


1. 


043 


c. 


4p^ 


-0.098 


0.561 


0.759 


f . 0.57 


-0.959 


0. 


361 


-c. 


If 8 


2.517 


-1.064 


n . 057 


1.626 


0.381 


-1. 


00 4 


0. 


886 


-0.063 


0.35 8 


-0.059 


0.381 


2.113 


-0 • 


518 


0. 


151 


-0. 150 


1.043 


0. 361 


-1.C04 


-0.518 


1. 


619 



T TEST OF THE HYPOTHESIS THAT A 
COMPONENT CF THE P2 VECTOR ECUALS 
7 F R C AT A .If SIGNIFICANCE LEVFL 



TF ST STATISTIC RESULT 

5.91295 REJECT 

6.03753 REJECT 

-C. 22025 ACCEPT 

1.44270 ACCEPT 

10.36179 REJECT 

-1.23196 ACCEPT 

2.79636 REJECT 
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ENGINE NUMBER: 524646 



DATA : 



Y 1 


Y2 


Y 3 


Y 4 


Y5 


Y 6 


Y7 


XI 


12. ' 


38.0 


2.0 


r . 1 


1 l. p 


2.0 


4.0 


32.0 


9. 


A 0 


7.r 


0 . 0 


13.0 


2.C 


6 • C 


32. n 


i c . : 


4 3 . 0 


1 . 0 


2. 0 


11.0 


2.0 


5.0 


27. 9 


9." 


47. 1 


0 . ^ 


0.0 


9.0 


3.0 


2.0 


22.0 


8.9 


3 0 . o 


2. U 


0.0 


10.C 


l.C 


l.G 


9.C 


7* ^ 


2 9.0 


c. : 


1. 7 


8.0 


1.0 


1.0 


8.0 


ll.o 


4 2.0 


3.0 


0.0 


12.0 


2 • 0 


5.0 


63. o 


11.0 


4 5.0 


6 • n 


o.n 


12.0 


2.o 


3 . C 


72.0 


15* r 


56.0 


• \ 


C.o 


19. n 


5 . 0 


7.C 


12 6.0 



NUMBER PE DATA FCINTS = 9 



FESLLTS : 



rSTIMCTt CF P 
*1 B2 



7.763C4 
31.59843 
2. 17532 
C.5P916 
F. 3 6977 
1 • 1 0 3 5 c 
2*06722 



''•05629 
G. 19645 
0.C9365 

^ 5 g n 

562 



’ *4 ft V 






C .02581 
0,03947 



ESTIMATE CF COVARIANCE MATRIX $IG V A 



1* 295 


1.34 3 


-3.71Q 


O.C 23 


2.37 7 


0. 26C 


C.69 3 


1. 345 


13.200 


C. 202 


0.146 


0.74 7 


1.910 


2. 558 


-C. 719 


0.202 


7.492 


-n,679 


3.54 6 


-1. 1C7 


1.287 


C.C23 


0 • 146 


-0,679 


9.514 


0 • C 2 0 


G . 0 1 3 


o. 336 


C. 377 


0.747 


C.546 


O.C 20 


2 • n 1 9 


°. 31 2 


1.564 


C.260 


1.910 


-1.107 


O.C 13 


0.312 


0. 554 


0. 244 


C.693 


2. 558 


1.287 


0.3 36 


1.564 


0.244 


2. 80 l 



T TEST CF THE HYPOTHESIS THAT A 
COMPONENT CF THE B2 VECTOR ECUALS 
7EPC AT A .1C SIGNIFICANCE LEVEL 



TFST STATISTIC 


RESULT 


5.30951 


REJECT 


5. 80300 


REJECT 


C . 1 4 1 1 0 


ACCEPT 


-o • 89370 


ACCEPT 


5.71138 


REJECT 


3. 72125 


REJECT 


2.53127 


REJECT 
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ENGINE NUMBER: 516627 



DA TA : 



Y 1 


Y 2 


Y3 


Y4 


Y5 


Y 6 


Y7 


XI 


e. ; 


32.0 


9.0 


0. 0 


7. r 


3.9 


3.0 


24. C 


4. 1 


18.9 


1.0 


0.0 


7. C 


3.0 


3.0 


10. o 


11.0 


37.9 


9.'' 


2.0 


14. C 


3.0 


4.0 


86.0 


9.0 


27.0 


2.0 


0.0 


9.0 


2 • C 


4.0 


24.0 


5.0 


2 3.0 


2.0 


1 . 9 


7.0 


1.0 


o.o 


3.0 


1C.C 


32.0 


3.0 


ri o 

'V • 0 


13.0 


2. r 


3.0 


6 3.0 


1 0 . r- 


39.0 


2.0 


0. 0 


13.0 


2.0 


A .o 


6 3. 0 


8.0 


3 5. 9 


o.o 


0.0 


13.0 


2.0 


3.0 


68. o 


6 . v 


23. 9 


2.° 


1.0 


11. o 


1.0 


l.o 


23.0 


6.0 


27.0 


2.0 


0.0 


ll.o 


1.0 


1.0 


36. r 


6.0 


30. * 


2. C 


0.9 


12.0 


3.0 


4.0 


60. r 



NUMBER PF PATA FCINTS = 11 



RESULTS: 



ESTIMATE cf p 



PI 

4 • 9547C 
21.32212 
2. ''96 7 5 
C . ? r 373 
6.62466 
1.672C2 
1.61738 



B2 

9.0612° 
0.1B915 
2.02352 
0. 00378 
0.09017 
0 . 0 0 ° 9 1 
0.01765 



2. 669 3.273 

3.27P 11.578 

2.382 3.910 

C . 1 1 6 0 . C 2 9 

0.339 -0.452 

0.021 -0.571 

0.172 -2.C98 



3.919 


0.029 


9.631 


o. 995 


C . 9 9 5 


0. 493 


-1.871 


0.094 


1.12 5 


-9.072 


1.274 


-0. 107 



-J.021 
-0.571 
1.125 
094 -0.072 

1.157 -0.718 

-0.718 0.684 

-0.746 0.922 



F ST I MATE CF CCVARIANCE MATRIX SIC-MA 



2.282 



C. 1 16 



-0.339 

-0.452 

-1.871 



(; 



C . 17 2 
-2. 088 
1. 274 
-C. 107 
-0. 746 
r. 922 
2.463 



T TEST OF THE HYPOTHESIS THAT A 
COMPONENT CF THE B2 VECTOR EQUALS 



ZERO AT A .1C 



TEST STATISTIC 

3.27932 
4. 85956 
C. 66246 
0.47103 
7. 32913 
1. 04718 
0.93349 



SIGNIFICANCE LEVEL 



RESULT 

REJECT 

REJFCT 

ACCEPT 

ACCEPT 

REJECT 

ACCEPT 

ACCEPT 
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ENGINE NIJMPFR: 516101 



TATA : 



Y 1 


Y 2 


Y3 


Y4 


Y5 


Y 6 


Y7 


XI 


16.*’ 


4 4. 0 


3. 0 


0.0 


23. C 


3.0 


4.C 


8 0. C 


16. T 


40. 0 


0.0 


0. 0 


23.0 


1.0 


l.C 


128.0 


5. ■- 


17.0 


l.n 


3.0 


8.0 


o.e 


4.9 


12 8.0 


1 r . 0 


33. 7 


9* G 


1. 0 


19.0 


1.0 


l.C 


192. ri 


1 1 . •'■ 


36.0 


3.0 


4 • / 


21. r 


1 ■ 0 


2.0 


114. 0 


6. . 


? 9 . 3 


9.C 


o.o 


17. f 


1 • 0 


f! . n 


56.0 


H .f 


3 7.0 


1.0 


C. J 


16.0 


2.C 


3.-: 


90.0 


7. » 


21.9 


5.R 


2.0 




1.0 


o . 0 


23. ° 


6. " 


2*t. 0 


2.0 


O.o 


ii.o 


1.0 


1.0 


23. C 


7.0 


2 2.0 


8.0 


2.0 


9.0 


l.C 


4 . C 


11.0 


11/ 


2 5 . j 


7.0 


o. 0 


1 5 . o 


1.0 


1.0 


35. o 



NUMBER r>p DATA PCI NTS = II 



PESUl TS : 



estimate TF R 
01 62 



P.f U12 
24.557 2‘ 
5. 14212 
r . 71736 
11.69666 
1. 276C 3 
1. 730 5 9 



0.C2141 
• U 6 5 7 6 
•0.0097 3 
o . one 12 

0.04Q 24 
•C.C01 18 
0.00223 



ESTIMATE CF CCVARIANCE MATRIX SIGMA 



13.487 
26.779 
-2. 4-6 c 
-2.767 
16. 654 
2. l. 4 ( 

r . 132 



26.779 
69. 565 
-6.069 
-7. 645 
38.864 



-2.465 
-6.069 
12.388 
0.348 
-0 . 676 



-2.767 
-7. 645 
0. 348 
1. 353 
-4. 921 



5.651 -0.344 -C .494 



i.676 -1.771 



16.654 
38. 864 
-0.676 
-4.921 
25.366 
2.733 



n » 6 3 5 -1.54 5 



2.04C 

5.651 

- 9.344 

-0.494 

2.733 

0.621 

0.363 



r . 1 3 2 

r. 676 

-1. 771 

r • 635 
-1. 545 
0. 36 3 
2. 75 C 



T TEST CF THE HYPOTHESIS THAT A 
CCMOONENT CF THE B2 VECTCP ECUALS 
7ERC AT A .10 SIGNIFICANCE LEVEL 



TEST STATISTIC 

1.04745 
1.41632 
-C. 4 9665 
( .0 1914 
1. 75634 
-0. 26838 
0.24171 



RESULT 

ACCEPT 
ACCEPT 
ACCEPT 
ACC E°T 
ACCEPT 
ACCEPT 
ACCEPT 
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ENGINE NUMBER 5 516132 



DATA : 



Y 1 


Y 2 


Y 3 


Y4 


7. " 


25.0 


9.0 


2.0 


16. ' 


29.0 


3.0 


O.o 


c. ^ 


35.0 


6.0 


2.o 


5. > 


31.0 


7 • C 


1.0 


5.' 


3 1 . 0 


4 . r 


0.0 


1 C. . 0 


2 9. 0 


7.0 


0.0 


6.i. 


23.0 


3.0 


0.0 


7.0 


24. 0 


2.0 


0.0 


7." 


2 '. 1 


10.0 


0. 0 


LUMBER OF 


C A T A 


POINTS = 


5 



Y5 


Y 6 


Y7 


XI 


B.r 


1.0 


2 • C 


26.0 


9.0 


2. 0 


a r 


100.0 


11.0 


2.0 


4.0 


139. 0 


9 . r 


2.0 


2.0 


7 3. 0 


9.0 


2.0 


C.C 


107. 9 


11.' 


1.0 


o . c 


5P.n 


7.0 


2 . C 


2.C 


44.0 


7.o 


1.0 


1.0 


58. 0 


5.0 


1 . 0 


2.0 


3 5.0. 



RESULTS 



ESTIMATE CF B 



PI 



R2 



5.56531 
19. 6876* 
7. 67337 
C . 253C e 
6 . C' 2 4 3 c 
C. 97049 
1.15398 



v. 04045 
0. 10 ROB 
-0.02335 
0.00425 
0.03403 
0 . 0 0 G 6 3 
0.004TB 



ESTIMATE CF CCVARIANC5 MATRIX SIGMA 



9.420 


-1.576 


-1.412 


-1.020 


0.551 


-0.225 


-2. 


660 


-1.576 


6.78? 


2.538 


1. 2 29 


3.546 


0.0 0 8 


C. 


036 


-1.412 


2.533 


7.866 


1.430 


0.868 


-0.327 


1. 


660 


-l.C 2C 


1.229 


1.4 30 


n • 86 r 


0.452 


-0.034 


0. 


941 


0.551 


3. 546 


0.86 8 


o.452 


2.472 


-C. 126 


-0. 


333 


-0.225 


0.CC8 


-C.327 


-0.034 


-0.126 


0. 170 


o. 


04 P 


-2.660 


0.036 


1.660 


0.941 


-0.333 


0.048 


2. 


00 5 



T TEST CF THE HYPOTHESIS THAT A 
COMPONENT CF THE 62 VECTOR EQUALS 
Z P R r AT A .10 SIGNIFICANCE LEVEL 



TEST STATISTIC RESULT 

1.33313 ACCEPT 

4.42127 REJECT 

-1.06572 ACCEPT 

C • * 8443 ACCEPT 

2.28583 REJECT 

2.47053 REJECT 

C. 30463 ACCEPT 
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F 



CCMFARISCN CF EST IN AT ED PARAMETERS AMONG ENGINES 



EST IMATEC 



SERIAL # 

516726 
515355 
516678 
5 1 5 C 5 6 
5C5C90 
524735 
5163C5 
5 16C94 
516545 
51675C 
516575 
516534 
505C8 l 
516592 
515531 
5 1 56 C 8 
516380 
515515 
516731 
515447 
516C46 
516C63 
509064 
524646 
516627 
5 1610 1 
516132 



SERIAL it 

516726 
515355 
516676 
515C58 
5C5C9C 
524735 
516305 
516C94 
516545 
5 16 75 C 
516575 
5 16534 
509081 
516552 
5 1593 1 
515608 
51638C 
515515 
516731 
515447 
516048 
5 1606 3 
505C64 
524646 
516627 
516101 
516132 



PARAMETERS 



MEAN 

3. 80C0 

4.2222 
5 . 09 CO 

4.5556 
3.1818 
2.5714 

3.2222 

5.5556 
4. 1000 
4.6667 

4.2222 
5.7778 
4 . 1 G CO 
3.1111 
3.2000 

4.6667 
5.09C9 

4.3333 

6.2222 
3.9000 
4.0000 

3.3333 
5.8889 

2.3333 
3.0909 
4.3636 

5.6667 



MEAN 

1 . 60 CO 
1.0000 
C. 7778 
0.5556 
0.2727 
1.5714 
0.4444 
0.6667 
1.0000 
0.6667 
0.5556 
1.1111 
0.2000 
0.7778 
C.40CC 
0.6667 
1.0000 
0.3333 
1.4444 
C.8CC0 
1.1111 
C • 4444 
1.1111 
0.3333 
C . 36 36 
C . 72 73 
0.5556 



FCR CFRCMIUM 



VARIANCE 

7.7333 
10.6544 
7.2500 
7.2778 
6.1636 
13.3406 
12.4444 
11.5278 
8 . 7667 
5 . 7 5 C 0 
10. 1544 
9 .1944 
10 • 5444 
4.8611 
6.8444 
10. 75CC 
5.49C5 

8. 7500 
1.4444 
5.8778 
6.5000 

ll.OCCC 

7.6111 

6.7500 
5.05C5 

11.4545 

8.CCC0 



FCR SILVER 



VARIANCE 

0.7111 

1.2500 
0.4444 
0.2778 
0.2182 
2.4176 
0.5278 
C. 75C0 
1.1111 
C. 5CC0 
C.5278 
1.6111 
0.1778 
0.4444 
0.2667 
0.5CC0 
1 .OCCO 
C.25CC 
0.7778 
1.0667 
0.6111 
0.2778 
C. 8611 
C.5CC0 
0.4545 
1.2182 
0.7778 



ESTIMATED PARAMETERS 
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ESTIMATED PARAMETERS FOR MAGNES. 



SERIAL H 


MEAN 


VARIANCE 


516726 


1. 6000 


0.4889 


515355 


1. 3333 


0.7500 


516678 


2.2222 


1. 1944 


515058 


0. 7778 


0.6944 


509090 


1.1818 


0.9636 


524739 


2.5714 


2. 2637 


516309 


0.6667 


1.7500 


516C94 


1.3333 


0.2500 


516949 


2. 1000 


0. 5444 


516750 


1.2222 


G. 9444 


516575 


2.7778 


3. 1944 


516534 


1.8889 


0. 3611 


509081 


1.6000 


1.3778 


5165 Q 2 


2. 3333 


0.7500 


515931 


1. 1000 


1.4333 


515608 


1.4444 


0.7778 


516380 


1. 0900 


G.C9C9 


515515 


2.0000 


0.2500 


516731 


1.4444 


0.7778 


515447 


2.3000 


3.7889 


516C48 


1.0000 
2. 3889 


0. 5000 


516063 


8.3611 


509064 


2. 6667 


2.2500 


524646 


2.2222 


1.4444 


516627 


2. 0^09 


0.6909 


516101 


1.1818 


0. 5636 


516132 


1.5556 


0.2778 


ESTIMATED 


PARAMETERS 


FOR NICKEL 


SERIAL n 


MEAN 


VAR I ANCE 


516726 


2.0000 


1.5556 


515355 


2.0000 


2.0009 


516678 


3.3333 


3.7509 


515058 


1.4444 


1.2778 


509090 


1.0909 


1.09C9 


524739 


2.9296 


2.9945 


516309 


1.0000 


1. 2509 


5 1 6094 


2. 3333 


3.5093 


516949 


3.20C0 


2.8444 


516750 


2. 5556 


2.2778 


516575 


2.8889 


3.3611 


516534 


3. 4444 


2.2778 


509C81 


0. 9000 


0 • 5444 


516592 


2. 6667 * 


3.2507 


515931 


2.3OC0 


1.5667 


515608 


3.1111 


3.6111 


516383 


2.0000 


3.4000 


515515 


3. 3333 


4. 0000 


516731 


3. 0000 


5.0000 


515447 


1.4000 


1.6000 


516048 


1.7778 


1.9444 


516063 


2.888° 


1. 1111 


509064 


4.0000 


3.0000 


524646 


3. 7778 


4.6944 


516627 


2. 3636 


2.4545 


5 1 6 1 C 1 


1. 9091 


2.4909 


516132 


1. 4444 


1.7778 
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REGRESSION LINES FOR ALUMINUM 



SERIAL # 


INTERCEPT 


SLOPE 


RESIDUAL 


SAMPLE 








VAR I ANCE 


CORRELATION 


516726 


8.7000 


0.61562 


2.5416 


0.6655 


515355 


7. 6667 


0 . losoe 


5.7870 


0.5887 


516678 


7 • 4444 


0. U205 


3.8633 


0.2^41 


515058 


7.4444 


-0. 31984 


17.5090 


-C. 1703 


509090 


14.0000 


0. 03755 


14.9953 


0.5162 


524739 


12.3571 


0.00456 


15.6695 


0.079C 


51 63C9 


6. 44.44 


0.07506 


6.4149 


0.7303 


5 16C C 4 


7.2222 


0. 03688 


4.6820 


0.4141 


5 1 6949 


6.7000 


-0.02339 


4.8Q07 


-0.3358 


516750 


6.1111 


C. 06514 


5.5361 


0.0591 


516575 


10. 5556 


0.03202 


27.8192 


0.1924 


516534 


14. 6667 


0.02293 


4. 1051 


0.3194 


5 C 90 8 1 


7 . COOC 


-0.00015 


25.4999 


-0.0019 


5 1 65 Q 2 


6.2222 


0.02678 


3 • 6684 


0.3622 


515931 


9.2000 


0.04747 


13.7138 


0.4366 


515608 


6.2222 


-0.00061 


7.0790 


-0.C074 


516380 


6.4845 


0.08528 


16.5446 


0.4085 


515515 


13.8889 


0. 08896 


11.6934 


0.5362 


516731 


3.4444 


-0.03564 


11.9624 


-0.3033 


5 1 5447 


3.4000 


0. 11910 


4.4045 


0. 74° 5 


516C48 


6.0010 


0.06004 


4.1227 


C • 731 1 


5 1 60 63 


5.5556 


0.03312 


6.3866 


C. 3314 


5C9064 


11.2222 


0.02695 


C.9903 


0.9128 


524646 


10.2222 


0.05629 


1. 2945 


0. 8950 


516627 


7.5455 


0.06129 


2.6692 


0.7378 


516101 


9.7273 


0.92141 


13.4872 


0.3296 


516132 


8.4444 


0. 04049 


0 • 4 20 0 


0.4659 




REGRESSION 


LINES FOR 


IRON 




SERIAL * 


INTERCEPT 


SLOPE 


RESIDUAL 


SAMPLE 








VAR I ANCE 


CORRELATICN 



516726 
515355 
516678 
515058 
5C90<?(‘ 
524739 
5163C5 
516054 
516944 
516750 
516575 
516534 
5C9C81 
516592 
515931 
515608 
51638C 
515515 
536731 
515447 
516048 
516C 63 
5C°C64 
524646 
516627 
5 1 6 1 Cl 
516132 



30. 1000 

27.3333 

26.5555 

22.6667 
32.4545 
33.6429 

26.4444 

26.6667 
24. 7000 

17.3333 

35.5555 

34.4444 
25. 50u0 

17.4444 
30.1000 

31.6667 
24. 3636 

39.3333 
21.1111 
29.2000 

34.2222 

19.2222 
35.7778 
40.1111 
29.8182 
29.8132 

27.4444 



0. 06875 
f. 11819 
0.09537 
0. 03226 
C. 10673 
0.06821 
C. 34263 
0.20229 
0. 17268 
0.13478 
0.19532 
0.18827 
0. 0587Q 
0.07297 
0. 11599 
0. 02774 
0. 17684 
0. 16124 
0.04852 
G. 22241 
0. 2245Q 
0. 12472 
0.07985 
0. 19645 
0.18915 
0.06576 
0.10908 



4.6265 
20.9136 
27.4270 
30.0517 
20.2778 
26.7694 
49.2843 
6.5651 
30.7814 
64.6504 
16.5041 
14. 60 30 
35.3305 
8.2043 
13. 5679 
39. 1967 
65.1 446 
20.4983 
17. 8794 
21.9142 
92.5444 
30.4751 
8.2291 
13. 1995 
11.5735 
69.5654 
6.7836 



0.9275 
0.7578 
0.0348 
0.2093 
0.3274 
0.6716 
0.8695 
0.9C34 
0.7165 
0.4137 
0.8408 
0.8263 
0.5359 
0.5779 
0.7662 
0.1417 
0.4236 
0.6562 
0.3342 
0.7159 
0.6459 
0.5130 
C. 9 169 
0.9099 
6. 8 50 9 
0.4269 
0.8581 
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REGRESSION LINES FOR COPPER 



SERIAL M 


INTERCEPT 


SLOPE 


RESIDUAL 

VARIANCE 


SAMPLE 

CORRELATION 


516726 


14. 5000 


0. 04810 


4.6190 


0.8667 


515355 


11.5556 


0.06015 


5.8640 


0.7449 


516678 


11.7778 


0.01226 


4.0435 


0.2029 


515058 


7.5556 


0.01855 


3.2891 


0.3^93 


*>r 9CQ0 


9.6364 


0.94157 


0.9376 


0.9364 


524739 


10.1429 


0.03660 


4. 3039 


0.7714 


516309 


8. 8889 


0.07843 


6.1235 


0.7527 


516094 


9.2222 


0.02886 


10.4862 


0.2314 


516949 


9. 6000 


0.06982 


2.4905 


0.8250 


516750 


10.7778 


-0.06335 


5.5536 


-0.5888 


516575 


10.0000 


0.09021 


2.9364 


0.8620 


516534 


11.4444 


0.10755 


6.0699 


0.7928 


5C9C81 


5.7000 


0. 01973 


4.6695 


0.5058 


516592 


8.8839 


0.063OQ 


4.1955 


0.6503 


515931 


10.1000 


0.04496 


1.7140 


C.7927 


515608 


15.7778 


0.07357 


11.1454 


0.5799 


516380 


9.45^5 


0.08577 


1C. 0621 


0.4999 


515515 


11.2222 


0.08609 


1.8030 


0.8428 


516731 


10.3333 


0.02090 


11.0115 


0.1910 


515447 


10.2000 


C. 06485 


3.2416 


0.6137 


516048 


9.7778 


0.05309 


10.5203 


0.5102 


516C63 


10.1111 


C. 04695 


6.8290 


0.4339 


509064 


14.0000 


0.06053 


1.6263 


0.9699 


524646 


11.6667 


0.07562 


2.0192 


0.9074 


516627 


10.6364 


0.09017 


1.1567 


0.9255 


5 161 Cl 


15. 6364 


0.04924 


25. 3663 


0.5052 


516132 


8 • 4444 


0.03403 


2.4722 


0.6539 
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